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Abstract 

We develop a new systematic procedure for the Regge limit in perturbative QCD to ar- 
bitrary logarithmic order. The formalism relies on the IR structure and the gauge symmetry 
of the theory. We identify leading regions in loop momentum space responsible for the sin- 
gular structure of the amplitudes and perform power counting to determine the strength of 
these divergences. Using a factorization procedure introduced by Sen, we derive a sum of 
convolutions in transverse momentum space over soft and jet functions, which approximate 
the amplitude up to power-suppressed corrections. A set of evolution equations generalizing 
the BFKL equation and controlling the high energy behavior of the amplitudes to arbitrary 
logarithmic accuracy is derived. The general method is illustrated in the case of leading 
logarithmic gluon reggeization and BFKL equation. 
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1 Introduction 



The study of semihard processes within the framework of gauge quantum field theories has a long 
history. For reviews see Refs. PP-P]- The defining feature of such processes is that they involve 
two or more hard scales, compared to Aqcd, which are strongly ordered relative to each other. 
The perturbative expansions of scattering amplitudes for these processes must be resummed since 
they contain logarithmic enhancements due to large ratios of the scales involved. One of the most 
important examples is elastic 2^2 particle scattering in the Regge limit, s 3> |t| (with s and t 
the usual Mandelstam variables). It is this process that we investigate in this paper. We extend 
the techniques developed in Refs. ^ and and devise a new systematic method for evaluation 
of QCD scattering amplitudes in the Regge limit to arbitrary logarithmic accuracy. 

The problem of the Regge limit in quantum field theory was first tackled in the case of fermion 
exchange amplitude within QED in Ref. 0. Here it was found that the positive signature ampli- 
tude takes a reggeized form up to the two loop level in Leading Logarithmic (LL) approximation. 
In Ref. the calculations were extended to higher loops, and the imaginary part of the Next- 
to-Leading Logarithms (NLL) was also obtained. The analysis in Refs. jU] and [7] was performed 
in Feynman gauge. It was realized in Ref. |H| that a suitable choice of gauge can simplify the 
class of diagrams contributing at LL. The common feature of all this work was the use of fixed 
order calculations. To verify that the pattern of low order calculations survives at higher orders, a 
method to demonstrate the Regge behavior of amplitudes to all orders is necessary. This analysis 
was provided by A. Sen in Ref. jl], in massive QED. Sen developed a systematic way to control 
the high energy behavior of fermion and photon exchange amplitudes to arbitrary logarithmic 
accuracy. The formalism relies heavily on the IR structure and gauge invariance of QED and 
provides a proof of the reggeization of a fermion at NLL to all orders in perturbation theory. 

The resummation of color singlet exchange amplitudes in non-abelian gauge theories in LL 
was achieved in the pioneering work of Ref. jU], where the reggeization of a gluon in LL was also 
demonstrated. The evolution equations resumming LL in the case of three gluon exchange was 
derived in Ref. [TO]. In Ref. fS], n-gluon exchange amplitudes in QCD at LL level were studied and 
a set of evolution equations governing the high energy behavior of these amplitudes was obtained 
at LL. A different approach was undertaken in Ref. JT]. Here n —>■ m amplitudes were studied in 
SU(2) Higgs model with spontaneous symmetry breaking. Starting with the tree level amplitudes, 
an iterative procedure was developed, which generates a minimal set of terms in the perturbative 
expansion that have to be taken into account in order to satisfy the unitarity requirement of the 
theory. See also Ref. jT2|. The extension of the BFKL formalism to NLL spanned over a decade. 
For a review see Ref. J^l- The building blocks of NLL BFKL are the emissions of two gluons or 
two quarks along the ladder, Ref. jHj, one loop corrections to the emission of a gluon along the 
ladder, Ref. [1^], and the two loop gluon trajectory, Refs. [12], [IZ], [TH] and [IH]. The particular 
results were put together in Ref. [20j. In Ref. the trajectory for the fermion at NLL was 

evaluated by taking the Regge limit of the explicit two loop partonic amplitudes, Ref. [22]- 

Besides the NLO perturbative corrections to the BFKL kernel a variety of approaches have 
been developed for unitarization corrections, Refs. [23112111235 which extend the BFKL formalism 
by incorporating selected higher-order corrections. The procedure proposed in this paper, in a way, 
places these approaches in an even more general context. In principle, it makes it possible to find 
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the scattering amplitudes to arbitrary logarithmic accuracy and to determine the evolution kernels 
to arbitrary fixed order in the coupling constant. The formalism contains all color structures and, 
of course, the construction of the amplitude to any given level requires the computation of the 
kernels and the solution of the relevant equations. 

The paper is organized as follows. In Sec. |21we discuss the kinematics of the partonic process 
under study and the gauge used. In Sec. El we identify the leading regions in internal momentum 
space, which produce logarithmic enhancements in the perturbation series. After identifying these 
regions, we perform power counting to verify that the singularity structure of individual diagrams 
is at worst logarithmic. The leading regions lead to a factorized form for the amplitude (First 
Factorized Form). It consists of soft and jet functions, convoluted over soft loop momenta, which 
can still produce logarithms of s/\t\. In Sec. E] we study the properties of the jet functions 
appearing in the factorization formula for the amplitude. We show how the soft gluons can 
be factored from the jet functions. In Sec. El we demonstrate how to express systematically 
the amplitude as a convolution in transverse momenta. In this form all the large logarithms 
are organized in jet functions and the soft transverse momenta integrals do not introduce any 
logarithms of s/\t\ (Second Factorized Form). We derive evolution equations that enable us to 
control the high energy behavior of the scattering amplitudes. In Sec. El we illustrate the general 
methods valid to all logarithmic accuracy in the case of LL and NLL in the amplitude and we 
examine the evolution equations at LL. Some technical details are discussed in appendices 1X1 - 
IeI The first appendix treats power counting for regions of integration space where internal loop 
momenta become much larger than the momentum transfer. In Appendix^we illustrate the origin 
of special vertices encountered in the paper. In Appendix O we show a systematic expansion for 
the amplitude leading to the first factorized form. In Appendix Owe list the Feynman rules used 
throughout the text. Finally, in Appendix |El we demonstrate the origin of extra soft momenta 
configurations (Glauber region) which need to be considered in the analysis of amplitudes in the 
Regge hmit. 

2 Kinematics and Gauge 

We analyze the amplitude for the elastic scattering of massless quarks 



within the framework of perturbative QCD in the kinematic region s ^ —t {Regge limit), where 
•s = {pa + Pb)^ and t = are the usual Mandelstam variables. We stress, however, that the 
results obtained below apply to arbitrary elastic two-to-two partonic process. We pick process (|H) 
for concreteness only. The arguments in Eq. label the momenta and the colors of the quarks 
(we do not exhibit the dependence on the polarizations). We choose to work in the center-of-mass 
(cm.) where the momenta of the incoming quarks and the momentum transfer have the following 
components ^ 



q{PA, ta) + q'{PB, Tb) q{pA - q, ri) + q'ips + g, ^2), 



(1) 



PA 




We use light-cone coordinates, v — (w+, v ,v±), v'- 



(^;0±^;3)/^/2. 
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9 = (0+.0-,gx). (2) 

Strictly speaking = ±\t\/\/2s, so the components vanish in the Regge hmit only. 
In the color basis 

with Nc the number of colors, we can view the amplitude for process ([[]) as a two dimensional 
vector in color space 



where Ai and As are defined by the expansion 

rs, ri r2 ~ ^1 (^l)rA rg, rir2 ~^ Ag (&8)r^ ^b, ri r2 • (5) 

Since the amplitude is dimensionless and all particles are massless, its components can depend, in 
general, on the following variables 

Ai = Ai ( Aj,^,as{iJ,'^),e) for i = 1, 8, (6) 

where /x is a scale introduced by regularization. We use dimensional regularization in order 
to regulate both infrared (IR) and ultraviolet (UV) divergences with D = 4 — 2e the number 
of dimensions. Choosing the scale fi^ = s, the strong coupling, is small. However, in 

general, an individual Feynman diagram contributing to the process at r-loop order can give 
a contribution as singular as (s/t) al'^Hn'^^{—s/t). In Sec. 15.31 we will confirm that there is 
a cancellation of all terms proportional to the i-th logarithmic power for i = r + 1, . . . ,2r at 
order a^"*"^ in the perturbative expansion of the amplitude. Hence at r loops the amplitude is 
enhanced by a factor (s/t) a^'^'^ lif{—s/t), at most. In order to get reliable results in perturbation 
theory we must, nevertheless, resum these large contributions. In the fc-th nonleading logarithmic 
approximation one needs to resum all the terms proportional to (s/t) a^"*"^ ln^~-^(— s/t), j = 
0, . . . , /c at r-loop level. 

We perform our analysis in the Coulomb gauge, where the propagator of a gluon with momen- 
tum k has the form 

Naf3{k,k) If kakfs + kakfs - kakp\ 
" -^^"^ W^e —k ) ' 

in terms of the vector 

k = k - {k-r])7], (8) 
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with 

an auxiliary four- vector defined in the partonic cm. frame. The numerator of the gluon propagator 
satisfies the following identities 

k ■ k 

k'^N^p{k,k) = 0. (10) 

The first equality in Eq. (fTUI) is the statement that the nonphysical degrees of freedom do not 
propagate in this gauge. For use below, we hst the components of the gluon propagator: 



N+-{k) = N-+{k) 
N++{k) = N~~{k) 



IXi IXi _L 

k+k- 



k-k 



N^\k)=N^^{k) = ±^^~ 



2k -k 

N'^(k) = N^'ik) = g'i (11) 

k ■ k 

We note that these are symmetric functions under the transformation k^ —k^, except for 
the components N^^ = N^^, which are antisymmetric under this transformation. It was demon- 
strated in Ref. |26] that QCD is renormalizable in Coulomb gauge, by considering a class of gauges 
which interpolates between the covariant (Landau) and the physical (Coulomb) gauge. 



3 Leading Regions, Power Counting 

In order to resum the Regge logarithms, we need to identify the regions of integration in the loop 
momentum space that give rise to singularities in the limit t/s 0. We follow the method devel- 
oped in Refs. [2Z1I2H1! which begins with the identification of the relevant regions in momentum 
space. 

3.1 Singular contributions and reduced diagrams 

The singular contributions of a Feynman integral come from the points in loop momentum space 
where the integrand becomes singular due to the vanishing of propagator denominators. However, 
in order to give a true singularity the integration variables must be trapped at such a singular 
point. Otherwise we can deform the integration contour away from the dangerous region. These 
singular points are called pinch singular points. They can be identified with the following regions 
of integration in momentum space, 

1. soft momenta, with scaling behavior A;^ ~ a ^/s for all components (cr ^ 1), 
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Figure 1: The reduced diagrams a) and c) contributing to the amphtude. Diagram b) represents 
a decomposition of diagram a) for the purpose of power counting. 



2. momenta colhnear to the momenta of the external particles, with scaling behavior 

~ ^/s, k~ X Y^s, ~ A^/^ for the particles moving in the + direction and 
A;''" ~ A ^/s, ~ ^/s, \k±\ A^^^ ^/s for the particles moving in the — direction, 

3. so-called Glauber or Coulomb momenta, Ref. [SH], with scaling behavior k^ ~ ^/s, 
\kj_\ ~ cr y/s, where A < cr^ < a, and where the scaling factors X,a satisfy the strong 
ordering A <^ a <^ 1 (The origin of this region is illustrated in Appendix lEl). 

4. hard momenta, having the scaling behavior k^^ ~ ^/s for all components. 

The extra gauge denominators l/{k ■ k) originating from the numerators of the gluon propagator, 
Eq. (jZj), do not alter the classification of the pinch singular points mentioned above. Actually, 
only the subsets 1 and 3 in the above classification can be produced due to the extra gauge 
denominators. 

With every pinch singular point, we may associate a reduced diagram, which is obtained from 
the original diagram by contracting all hard lines (subset 4) at the particular singular point. As 
shown in Refs. [23l2HlEni the reduced diagram corresponding to a given pinch singular point must 
describe a real physical process, with each vertex of the reduced diagram representing a real space- 
time point. This physical interpretation suggests two types of reduced diagrams contributing to 
the process (P), shown in Fig. [H 

The jet A{B) contains lines whose momenta represent motion in the + (— ) direction. The 
lines included in the blob 5" and the lines coming out of it are all soft (configurations 1 and 3 in 
the classification of loop momenta described above). These two oppositely moving (virtual) jets 
may interact through the exchange of soft lines. Fig. and/or they can meet at one or more 
space-time points. Fig. 

Having found the most general reduced diagrams giving the leading behavior of the amplitude 
for process (Q) in the Regge limit, we can estimate the strength of the IR divergence of the integral 
near a given pinch singular point. First we restrict ourselves to cases involving subsets 1 and 2 
from the classification of loop momenta above. To do so, we count powers in the scaling variables 
A and a. 
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The scaling behavior of these loop momenta implies that every soft loop momentum contributes 
a factor cx^, every jet loop momentum gives rise to the power A^, every internal soft boson (fermion) 
line provides a contribution cx"^ and every internal jet line (fermionic or bosonic) scales as 

A~^. In addition, there can be suppression factors arising from the numerators of the propagators 
associated with internal lines and from internal vertices. As pointed out in Ref. |2ZI; in physical 
gauges each three-point vertex connecting three jet lines is associated with a numerator factor 
that vanishes at least linearly in the components of the transverse jet momenta, and therefore 
provides a suppression A^''^. 

We are now ready to estimate the power of divergence corresponding to the reduced diagrams 
describing our process. First we restrict ourselves to the case shown in Fig. [T^. As indicated 
schematically in Fig. ^p, we can perform the power counting for the jets and for the soft part 
separately. All soft propagators and all soft loop momenta are included in the soft subdiagram S. 
The superficial degree of IR divergence of the reduced diagram R from Fig. QJi and Fig. can 
then be written as 

u{R) = iu{A) + iu{B) + iu{S), (12) 

where the external lines and loops of S' are included in S. For uj{R) > the overall integral is 
finite, while ll!{R) < corresponds to an IR divergent integral. When uj{R) = 0, the integral 
diverges logarithmically. Here we set A ~ cr for power counting purposes. We come back to the 
effect of relaxing this condition in connection with a discussion of item 3, Glauber regions, in our 
list of singular momentum configurations. 



3.2 Power counting 

In this subsection, we consider the case when all vertices in a diagram are elementary only, that 
is, without contracted subdiagrams carrying large loop momenta. In Appendix^ we show that 
our conclusions are unchanged by contracted vertices. 

We perform the power counting for the soft part S first. Let f,b be the number of fermion, 
boson lines external to S and let E = f + b. The superficial degree of divergence for S, found by 
summing powers of a, can be written 

uj{S) =A{E -2) -2b~ f + 2 + uj{S'), (13) 

where the first term is due to loop integrations linking S' to the jets, while the second and the third 
terms originate from propagators associated with the bosonic and fermionic lines, respectively, 
connecting the jets A, B and the soft part S' . The term +2 is introduced because we are resumming 
only leading power corrections proportional to s/t and therefore we exclude the overall factor s/t 
from the power counting. Since the lines entering 5" are soft, we obtain the superficial degree of 
divergence for 5" simply from dimensional analysis. It is given by 

oj{S')=A-b-?,f/2. (14) 

Combining Eqs. (fT^ and (jHj), the superficial degree of infrared divergence for the soft part S is 
then 

cj(5) =6 + 3//2-2. (15) 
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Before carrying out the jet power counting, we introduce some notation. Let Ea be the number 
of soft hues attached to jet A; I is the total number of jet internal lines; is the number of a- 
point vertices connecting jet lines only; Wa has a meaning similar to Va, with the difference that 
every vertex counted by Wa has at least one soft line attached to it. These are the vertices that 
connect the jet A to the soft part 5*. Finally, L denotes the number of loops internal to jet 
A. As noted above, we will perform the power counting for the case when the scaling factor for 
the soft momenta, a, is of the same order as the scaling factor for jet A momenta. When the 
scaling factors are different we encounter subdivergencies, which can be analyzed the same way 
as described below. We also assume that there are no internal and external ghost lines included 
in the jet function. Later we will discuss the effect of adding ghost lines. 

The superficial degree of divergence for jet A can now be expressed as 

u{A) = 2L- I + V3/2. (16) 

The last term represents the suppression factor associated with the three point vertices. We denote 
the total number of vertices internal to jet A by 

V = '^iVa+Wa). (17) 

a 

Next we use the Euler identity relating the number of loops, internal lines and vertices of jet A 

L = I-v + l, (18) 
and the relation between the number of lines and the number of vertices 

21 + EA + 2 = J2c^iva + w^). (19) 

a 

Using Eqs. (jl6j) - (jl9j) we arrive at the following form for the superficial degree of divergence 
for jet y4 

u{A) = 1-{Ea + W3)/2 + ^(« - A){vo, + Wa)/2. (20) 

«>5 

Since every vertex counted by Wa connects at least one external soft line, we have the condition 

EA>W3 + J2wa. (21) 
a>4 

The equality holds when there is no vertex with two or more soft lines attached to it. Combining 
Eqs. ()20| ) -(PT |) we arrive at the following lower bound on the superficial degree of divergence for 
jet A: 

u{A) > l-EA + 5]]«^a/2 + 5^(a-4)(t;„ + w„)/2. (22) 

The third and the last term in Eq. (j22p are always positive or zero and hence 

iu{A) >1-Ea. (23) 



8 



A similar result holds for jet 5, and therefore the superficial degree of collinear divergence for jets 
A and B is 

uj{A) + u{B)>2-E, (24) 

with E = Ea + Eb as in Eq. (fT^ . Combining the results for soft and jet power counting, Eqs. 
ffT3j) and (I21|), respectively in Eq. (O, we finally obtain the superficial degree of IR divergence 
for the reduced diagram in Fig. [T^, 

u{R) > f/2. (25) 

This condition says that we can have at worst logarithmic divergences, provided no soft fermion 
lines are exchanged between the jets A and B. We can therefore conclude that a reduced dia- 
gram from Fig. ^ containing elementary vertices can give at worst logarithmic enhancements in 
perturbation theory. In order for the divergence to occur, the following set of conditions must be 
satisfied: 

1. There is an exchange of soft gluons between the jets A and B only, with no soft fermion 
lines attached to the jets. 

2. The jets A and B contain 3 and 4 point vertices only, see Eq. ()22|) . 

3. Soft gluons are connected to jets only through 3 point vertices, Eq. (f2^ . and at most one 
soft line is attached to each vertex inside the jets, Eq. (j^H). 

4. In the reasoning above we have assumed that there is no suppression factor associated with 
the vertices where soft and jet lines meet. In order for this to be true, the soft gluons must 
be connected to the jet A{B) lines via the +(— ) components of the vertices. 

Next we consider adding ghost lines to the jet functions. As we review in Appendix El the 
propagator for a ghost line with momentum k is proportional to l/{k ■ k). Hence every internal 
ghost line belonging to the jet gives a contribution which is power suppressed as 1/s. Since the 
numerator factors do not compensate for this suppression, we can immediately conclude that the 
jet functions cannot contain internal or external ghost lines at leading power. 

So far we have not taken into account the possibility when the soft loop momenta are pinched 
by the singularities of the jet lines. This situation allows different components of soft momenta 
to scale differently. For example, a minus component of soft momentum can scale as the minus 
component of jet A momentum A, while the rest of the soft momentum components may scale as 
a, where A -C a <C 1. The origin of these extra pinches is illustrated in Appendix lEI 

Let us see what happens when we attach the ends of a gluon line with this extra pinch to jet 
A at one end and the soft subdiagram S at the other end. The integration volume for this soft 
loop momentum scales as Ao"^. The soft gluon denominator gives a factor a~^. If this soft gluon 
is connected to the soft part at a 4-point vertex, there is no new denominator in the soft part. On 
the other hand, if the soft gluon is attached to the soft part via a 3-point vertex then the extra 
denominator including the numerator suppression factors scales as a~^. The new jet line scales 
as A"^ as long as the condition A^/^ ^ c" is obeyed; otherwise, we have the scaling u"^ for the 
extra jet line. For A^/^ > a the Glauber region produces logarithmic infrared divergence. When 
A^/^ < cr, the overall scaling factor A/cr^ indicates power suppressed contribution. 
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a) 



b) 



Figure 2: Jet A moving in the + direction (a) and jet B moving in the — direction (b). 

Let us now investigate another possibihty, when the soft gluon connects jet A and jet B directly 
and its momentum is pinched by the singularities of the jet A and the jet B lines. Denoting the 
scaling factors of jet A and jet B as Xa and A^, respectively, the integration volume provides the 
factor AaA^ct^ and the soft gluon denominator contributes the power cr^^. The extra jet A and jet 
B denominators scale as A^^ and A^^, provided A^^ > a and X^^"^ > a. For A^^ < a both extra 

jet denominators provide the scaling factor o"~^. When A^^ > cr, the power counting suggests 
logarithmically divergent integrals. 

We have therefore verified that when the softest component of a soft line satisfies the ordering 

~ ^ '^5 the Glauber (Coulomb) momenta produce logarithmically IR divergent integrals 
and need to be taken into an account when identifying enhancements in perturbation series. The 
analysis demonstrated above for the case of one Glauber gluon can be extended to the situation 
with arbitrary number of Glauber gluons. This follows from dimensional analysis, in a similar 
fashion as the treatment of purely soft loop momenta above. 

We conclude that the reduced diagram in Fig. ^ is at most logarithmically IR divergent, 
modulo the factor s/\t\. The reduced diagram in Fig. looses one small denominator compared 
to the reduced diagram in Fig. ^ and since we are working in physical gauge, this loss cannot 
be compensated by a large kinematical factor coming from the numerator. Hence the reduced 
diagram in Fig. [Hd is power suppressed compared to the reduced diagram in Fig. [T^, and we do 
not need to consider it at leading power. 

Finally, let us discuss the scale of the soft momenta. In the case of soft exchange lines, each 
gluon propagator supplies a factor l/(cr^ s), which we want to keep at or below the order t in the 
leading power approximation. Thus the size of the scale is fixed to be cr ~ \/s. In the case 
of soft lines which are attached to jet A or to jet B only, the scaling factor lies in the interval 
{\/\t\/y/s, 1). In the case of Glauber momenta, we again need a ~ \/R[/v^- Then the condition 
A^/^ > cr, which is necessary for the logarithmic enhancement, implies that the scaling factors for 
+ and — components of the Glauber (Coulomb) momenta can go down to \t\/s, the scale of the 
small components of jet momenta. Additionally, we should note that soft and jet subdiagrams 
that do not carry the momentum transfer may approach the mass shell (A, cr — *• 0). Such lines 
produce true infrared divergences, which we assume are made finite by dimensional regularization 
to preserve the gauge properties that we will use below. The same power counting as above shows 
that these divergences are also at worst logarithmic. 
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3.3 First factorized form 



The analysis of the previous subsection suggests the following decomposition of the leading reduced 
diagram from Fig. Let us denote the {n + 2)-point and (m + 2)-point Green functions, IPI in 
external soft gluon lines, corresponding to jet A, ■^(■^j^^ ^"^(pa, q, Vi^i, ■ ■ ■ , kn), Fig. |2h, and to jet 

J{B)ll'.'.'.lZ^PB,q,Tl]Pu ■ ■ ■ ,Pm), Fig. inb, respectively. The jet function J^^^j (JIb)) ^^^^ depends 
on the color of the incoming and outgoing partons ta, ti {tb, T2), as well as on their polarizations 
\a, Ai (Ab, A2), respectively. In order to avoid making the notation even more cumbersome we 
do not exhibit this dependence explicitly. In addition the dependence of J^^^j and jj^^^ on the 
renormalization scale n and the running coupling ^^(/i) is understood. The jet functions also 
depend on the following parameters: the gauge fixing vector 77, Eq. Q, of the Coulomb gauge, 
the four momenta of the external soft gluons attached to jet A (B), fci, . . . , fc„ (pi, • • • ,Pm), and 
the Lorentz and color indices of the soft gluons attached to the jet A (B), /xi, . . . , ai, . . . , a„ 
(i/i, . . . , z/^; bi, . . . , bm)- The momenta of the soft gluons attached to the jets A and B satisfy the 
constraints ^"^^ ki = q and Yl^=iPj — 9- 

According to the results of the power counting, the soft gluons couple to jet A via the minus 
components of their polarizations, and to jet B via the plus components of their polarizations. 
Therefore, only the following components survive in the leading power approximation 

j{m)bl...bm f ^ ^ „, . r, 

'Jb [PB,q,V,VA,Pl, 

where we have defined light-like momenta in the plus direction va = (1, 0, 0^) and in the minus 
direction vb = (0, 1,0^). We can now write the contribution to the reduced diagram in Fig. [T^, 
and hence to the amplitude for process (^J), in the form 

X S^^Sl^M-bj9,V,'VA,VB; ki,..., kn;pi, ...,Pm) 

X J^j^''''-'-{pB,q,V,VA;Pl,...,Pm), (27) 

where the sum over repeated color indices is understood. Corrections to Eq. ^I7\ are suppressed 
by positive powers of t/s. The jet functions Ja,b are defined in Eq. in the leading power 

accuracy. The internal loop momenta of the jets A, B and of the soft function S are integrated 
over. The soft function will, in general, include delta functions setting some of the momenta 
ki, . . . ,kn and color indices ai, . . . , of jet function Ja to the momenta pi, ■ ■ ■ ,Pm and to the 
color indices 61, . . . , &m of jet function Jb- The construction of the soft function S is described in 
Appendix O For a given Feynman diagram there exist many reduced diagrams of the type shown 
in Fig. and one has to be careful in systematically expanding this diagram into the terms 



• ; kn) 



,Pr, 



i=l / 

m \ 

n^^' (26) 
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that have the form of Eq. (|27|). This systematic method can be achieved using the "tuhp-garden" 
formahsm first introduced in Ref. [221 and used in a similar context in Ref. For convenience 
of the reader we summarize this procedure in Appendix O 

Let us now identify the potential sources of the enhancements in ln(s/|t|) of the amplitude 
given by Eq. (jTfj) . If we integrate over the internal momenta of the jet functions then we can 
get ln((p^ ■ /\t\) from Ja and ln((pB ■ /\t\) from Jb- In addition, according to the results of 
the power-counting, Eq. (j23j) . we know that the jet function with n external soft gluons diverges 
as 1/A"~^. After performing the integrals over the minus components of the external soft gluon 
lines attached to jet A and over the plus components of the external soft gluons connected to 
jet 5, these divergent factors are potentially converted into logarithms of ln((p^ ■ /\t\) and 
ln((p5 ■ /\t\), respectively. Our goal will be to separate the full amplitude into a convolution 
over parameters that do not introduce any further logarithms of the form ln(s/|t|). This task 
will be achieved in Sec. 15.11 In the following section, we analyze the characteristics of the jet 
functions. 



4 The Jet Functions 

In this section we study the properties of the jet functions A, B given by Eq. (PBj) since, as Eq. 

suggests, they will play an essential role in later analysis. Since the methods for both jet 
functions are similar we restrict our analysis to jet A only; jet B can be worked out in the same 
way. In Sec. 14.11 we examine the properties of jet A when the minus component of one of its 
external soft gluon momenta is of order \/\t\. In Sec. 14.21 we find the variation of jet A with 
respect to the gauge fixing vector r/, and finally in Sec. 14.31 we examine the dependence of jet A 
on the plus component of a soft gluon momentum attached to this jet. 



4.1 Decoupling of a soft gluon from a jet 

According to the results of power counting above, soft gluons attach to lines in jet A via the 
minus components of their polarization. Following the technique of Grammer and Yennie [SS! 
we decompose the vertex at which the jth gluon is connected to jet A. We start with a trivial 
rewriting of J a in Eq. ()26|) 




We now decompose the metric tensor into the form g^^ = K^^{kj) + G^'^{kj) where for a gluon 
with momentum kj attached to jet A, K^^ and G^^ are defined by 



va ■ kj — ie 

G^^'ikj) = g^" ~ K^'^'ikj). (29) 

The K gluon carries scalar polarization. Since the jet A function has no internal tulip-garden 
subtractions (they are contained in the soft function S), we can use the Ward identities of the 
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a) 



b) 



Figure 3: a) Decoupling of a K gluon from jet A. b) Leading contributions resulting from the 
attachment of a G gluon to jet A. 

theory [SI], which are readily derived from its underlying BRS symmetry to decouple this 
gluon from the rest of the jet A after we sum over all possible insertions of the gluon. The result 
is 

1 " 

jin)a,...a,...a„^^^^ Vb, T]' . . . , . . . , kj , . . . , = V {-igsP"''"') 

■ f^j — T-^ ... 

^ j{n l)ai...c,...a^. -''"(^p^^ T]] ki, . . . , ki + kj, ...,k-,... k^) ■ 

(30) 

The notation and k- indicates that the jet function J^""^"* does not depend on the color index 
and the momentum kj, because they have been factored out. In Eq. (jSU)) . Qs is the QCD coupling 
constant and f^^^-^"-! are the structure constants of the SU (3) algebra. The pictorial representation 
of this equation is shown in Fig. The arrow represents a scalar polarization and the double line 
stands for the eikonal line. The Feynman rules for the special vertices and the eikonal lines in Fig. 

are listed in Appendix O Strictly speaking the right-hand side of Eq. (j3Up and Fig. |2K contain 
contributions involving external ghost lines. However, from the power counting arguments of Sec. 
13.21 we know that when all lines inside of the jet are jet-like, the jet function can contain neither 
external nor internal ghost lines. Therefore Eq. (jHUIl is valid up to power suppressed corrections 
for this momentum configuration. 

The idea behind the K-G decomposition is that the contribution of the soft G gluon attached to 
the jet line in the leading power is proportional to v'^G^uVa = 0- order to avoid this suppression, 
the G gluon must be attached to a soft line. The general reduced diagram corresponding to the 
G gluon attached to jet A is depicted in Fig. Eb- The lines coming out of S as well as the lines 
included in it are soft. The letter G next to the jth gluon in Fig. reminds us that this gluon 
is a G-gluon attaching to jet J{A}^i via the G^^{kj) vertex. 

The reasoning described above applies to the case when all components of soft momenta are 
of the same order. In the situation of Coulomb (Glauber) momenta, this picture is not valid 
anymore, since the large ratio k±/k~ coming from the G"*"-*- component can compensate for the 
suppression due to the attachment of the G part to a jet A line via the transverse components of 
the vertex. 
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a) b) 
Figure 4: The result of a variation of jet function J^J^'' with respect to a gauge fixing vector. 

4.2 Variation of a jet function with respect to a gauge fixing vector 77 

In this subsection we find the variation of the jet function jj^"^ with respect to a gauge fixing 
vector 7]. The motivation to do this can be easily understood. We consider the jet function with 
one soft gluon attached to it only, J^^\pA,q,VB,v)- Let us define 

^a=Pa-V and (b = V '^b- (31) 

In these terms, jet function J^^^ can depend on the following kinematical combinations: 

Ja\pa, g, vb,v) = Ja\^a, Pa - vb, Cb, t). Using the identity Pa- vb = 2^aCb and the fact, that 
the dependence of Ja on the vector vb is introduced trivially via Eq. (|^. we conclude that 

j'A\pA,q,VB,v)=CBJ^A\U t). (32) 

Our aim is to resum the large logarithms of ln{p\) that appear in the perturbative expansion of 
the jet A function. In order to do so, we shall derive an evolution equation for p\ dJ^^ /dp\. Since 
Pa appears in combination with rj only, we can trace out the p\ dependence of J^^ by tracing out 
its dependence on rj. This can be achieved by varying the gauge fixing vector t]. The idea goes 
back to Collins and Soper [221 and Sen We will generalize the result to J^'' in Sec. 15.21 

We consider a variation that corresponds to an infinitesimal Lorentz boost in a positive + 
direction with velocity 6p. Thus, for the gauge fixing vector t] = (1,0,0,0) ^, Eq. (jHl), the 
variation is: 61] = fj 6f3 = (0, 0, 0, 1) 6f3. It leaves invariant the norm r^^ = 1 to order 0{6f3). The 
precise relation between the variation of the jet A function with respect to p\ and 61]°' is 

F) 7(1) f) 7(1) p, 7(1) f) 7(1) 

+ OJa_ _ -g O^A I ^ OJa _ ~a (y-JA , t(1) /ooN 

We have used the chain rule in the first equality and the simple relation (b dJ^^ /d(B = Ja\ 
following from Eq. ()H2j) . in the second one. 

In order for Eq. (jH^ to be useful, we need to know what the variation of jet A with respect 
to the gauge fixing vector rj is. The result of this variation for is shown in Fig. |3] It can be 

^For the moment we use Cartesian coordinates. 
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derived using either the formahsm of the effective action, Ref. or a diagrammatic approach 
first suggested in Ref. [S2] and performed in axial gauge. We give an argument how Fig. E] arises 
in Appendix El Here we only note that the form of the diagrams in Fig. 0] is a direct consequence 
of a IPI nature of the jet functions. The explicit form of the boxed vertex 

-iS''{k) = -i{r]-kfj'' + fi-k7]''), (34) 

as well as of the circled vertex is given in Fig. ^Jof Appendix |Dl while their origin is demonstrated 
in Appendix El The dashed lines in Fig. H] represent ghosts, and these are also given in Fig. 
of Appendix O The four vectors rj, given in Eq. ©, and 

appearing in Eq. (jH^ are defined in the partonic cm. frame, Eq. (jSJ. We list the components of 

S^N^'^ik) 

S,{k)mHk) = 

S,{k)m\k) = ^^k\ (36) 

for later reference. 

In Fig. 01 we sum over all external gluons. This is indicated by the sum over i. In addition, 
we sum over all possible insertions of external soft gluons {ii, . . . ,in^} G {!,... ,n}\{i}. This 
summation is denoted by the symbol vr. We note that at lowest order, with only a gluon i 
attached to the vertical blob in Fig. this vertical blob denotes the transverse tensor structure 
depending on the momentum ki of this gluon 

z (k^g-^ - ktk^) . (37) 

It is labeled by a gluon line which is crossed by two vertical lines. Fig. ^| The ghost line 
connecting the boxed and the circled vertices in Fig. can interact with jet A via the exchange 
of an arbitrary number of soft gluons. We do not show this possibility in Fig. ^ for brevity. 

Let us now examine what the important integration regions for a loop with momentum k in 
Fig. are. The presence of the ghost line and of the nonlocal boxed vertex requires that in the 
leading power the loop momentum k must be soft. It can be neither collinear nor hard. This will 
enable us to factor the gluon with momentum k from the rest of the jet according to the procedure 
described in Sec. 14.11 



4.3 Dependence of a jet function on the plus component of a soft 
gluon's momentum attached to it 

In this subsection we want to find the leading regions of the object kj'dJ^^^ /dkj' . This information 
will be essential for the analysis pursued in the next sections. For a given diagram contributing 
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a) 



b) 



Figure 5: a) Momentum flow of the external soft gluon inside of jet A. b) Typical contribution 
to k+d/^^/dk+. 

to we can always label the internal loop momenta in such a way that the momentum kj flows 
along a continuous path connecting the vertices where the momentum kj enters and leaves the 
jet function J^""*. When we apply the operation k^d/dk^ on a particular graph corresponding 

to Ja \ it only acts on the lines and vertices which form this path. The idea is illustrated in 
Fig. EK- The gluon with momentum k attaches to jet A via the three-point vertex vi. Then 
the momentum k flows through the path containing the vertices fi,f2,f3 and the lines /i,/2- 
The action of the operator k~^d/dk~^ on a line or vertex which carries jet-like momentum gives 
a negligible contribution, since the + component of this lines momentum will be insensitive to 
k~^. In order to get a non-negligible contribution, the corresponding line must be soft. In Fig. Et, 
lines h and I2 must be soft in order to get a non-suppressed contribution from the diagram after 
we apply the k~^d/dk~^ operation on it. This, with the fact that the external soft gluons carry 
soft momenta, also implies that the lines I3, . . . , Iq must be soft. This reasoning suggests that in 
general a typical contribution to kj'dJ^^ /dk^ comes from the conflgurations shown in Fig. 
It can be represented as 

J'a^"'-"" = /(^ifd^^^j j^-'-'^^^-"'<--<'ivA,q,V;h,...,k^;k[,...,k'^,) 

X /f'''^-'''-'{pA,q,V,VB;k[,...,k'^,). (38) 

The function contains the contributions from the soft part S and from the gluons connecting 

the jet ^ and S in Fig. Et". The jet function ■* has fewer loops than the original jet function 
Now applying the operation k^d/dk'^ to Eq. (jSHI), the operator kj'd/dkj' acts only to the 
function Hence we can write 

k^^jf^'-- = I (jl^'^k^kl^j^-'-'^^^--'^-^^^ 

j \ i=l / 3 

X jf^"'--<'(p^,g,r/,t;5;A;i,...,fc;,). (39) 
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We conclude that the contribution to kj~dJ^^ /dkj' can be expressed in terms of jet functions j'^ ■* 
which have fewer loops than the original jet function. 

5 Factorization and Evolution Equations 

We are now ready to obtain evolution equations which will enable us to resum the large logarithms. 
First, in Sec. 15.11 we will put Eq. ^I7\\ into what we call the second factorized form. Then, in Sec. 
15. 2[ we derive the desired evolution equations. In Sec. 15. 3t we will show the cancellation of the 
double logarithms and finally in Sec. 15. 4| we demonstrate that the evolution equations derived in 
Sec. 15.21 are sufficient to determine the high-energy behavior of the scattering amplitude. 

5.1 Second factorized form 

The goal of this subsection is to rewrite Eq. (j27jl into the following form ^ 

1=1 J \j=i J 

X Sj^'^l V, Va, Vb] . . . , kn±]PlL, Pm±] M) 

X TP''-''-{pB,q,7],VA;Pi±,...,Pm±;M), (40) 

where and F^^ are defined as the integrals of the jet functions and over the minus 
and plus components, respectively, of their external soft momenta, with the remaining light-cone 
components of soft momenta set to zero. 



n-l / pM \ 



1=1 
+ 



ki — 0, . . . , k^ — 0, k^ , . . . , k^ ), 

^^B^^^"'^"'{PB,q,V,VA]Pl±,---,Pm±;M) = YlyJ dpfj/j^^^^'"^"'{pB,q,V,VA]Pl±,---,PmL, 

pr = o,...,p;; = o,p^,...,p+). (4i) 

In Eq. (001) 5 S' is a calculable function of its arguments and M is an arbitrary scale of the order 
^/\t\. The functions Ta,b and S' depend individually on this scale, but the final result, of course, 
does not. Based on the discussion at the end of Sec. 13.31 one can immediately recognize that all the 
large logarithms are now contained in the functions F^ and F^. The convolution of F^, F^ and S' 
is over the transverse momenta of the exchanged soft gluons. Since these momenta are restricted 
to be of the order \/\t\., the integration over transverse momenta cannot introduce ln(s/|t|). This 
indicates that at leading logarithm approximation the factorized diagram with the exchange of 
one gluon only contributes. In general, when we consider a contribution to the amplitude at 
L = La + Lb + Ls' loop level, where La, Lb and Ls' is the number of loops in Ta,Tb and S', 
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respectively, we can get L — Ls' logarithms of s/\t\ at most. Hence, the investigation of the s/t 
dependence of the full amplitude reduces to the study of the and dependence of Ta and 
r^, respectively. We formalize this statement at the end of Sec. 15.31 after we have proved that Ta 
(Tb) contains one logarithm of (p^) per loop. 

Let us now show how we can systematically go from Eq. (jTTj) to Eq. (jlU)). We follow the 
method developed in Ref. We start from Eq. ^I7\ and consider the k~ integrals over the jet 
function J a for fixed k^, kij_: 

/n— 1 
n dkr R^^-^-{k^, ...,k-;...) J'^a^"'-"- {pa. g, V, vb; ^i, • • • , ^n), (42) 

7t i = l 

where Ra is given by the soft function S and the jet function Jb, 

m \j=l J 

X jt'^''^"'^"{PB,q,ri,VA\Pi,...,Pm)- (43) 
We next use the following identity for Ra- ^ 

n— 1 

i?^(A;r,...,A;-_i) = RA^k^ = . . . .k-_, = ^) \[e{M - \kr\) 

1=1 

n-l 



1) ■ ■ ■ ) 



kni Pli ■ ■ ■ 1 Pm) 



+ [RA{ki ,...,ki , fci+i = 0, . . . , = 0) 

i=l 

- RAik^, . . . , k7_^, A;r = 0, . . . , k'., = 0) ^(M - \kr\)] 

n-l 

X l[e{M-\kj\). (44) 
j=i+i 

We have suppressed the dependence on the color indices and other possible arguments in Ra for 
brevity. The scale M can be arbitrary, but, as above, we take it to be of the order of ^/\t\- The 
first term on the right hand side of Eq. (jl^ has all k~ = 0. The rest of the terms can be analyzed 
using the K-G decomposition discussed in Sec. 14.11 Consider the (z = 1) term, say, in the square 
bracket of Eq. ()44|) inserted in Eq. ()42|) . Let us denote it Ai. In the region | <C M the 
integrand vanishes. On the other hand, for \ki\ ~ M we can use the K-G decomposition for the 
gluon with momentum ki. The contribution from the K part factorizes and the integral over the 
component fcf has the form 

/dk~ 
[RT-'^'iK, k^ = 0,..., k-_, = 0) - 9iM - \k^\) Rl^-'^-ik^ = 0, . . . , = 
Va ■ fCl 

n-l / „Af \ 

X Y. I 'dsn'"'''' j J^r'^^'-^'-^^iPA^ q, V, vb; k2,...,k, + h,...,kn)j. (45) 

i=2 \ J— 2 / 

■^Recall that kn — q — {ki + . . . + kn-i), so kn is not an independent momentum. 
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Eq. (|45|) is valid when all the lines inside the jet are jet-like. In that case the contributions from 
the ghosts are power suppressed. The contribution corresponding to a G gluon comes from the 
region of integration shown in Fig. Eb- It can be expressed in the form of Eq. involving some 
j'j^ ■* with fewer loops than in the original J^^\ and an i?^ with more loops than in the original 
Ra- Then we can repeat the steps described above with this new integral. 

Every subsequent term in the square bracket of Eq. can be treated the same way as the 
first term. This allows us to express the integral in Eq. ()42|) in terms of k~ integrals over some 



J4" ■'s, which have the same or fewer number of loops than the original J^"^ 



r 



(n')ai.-<, 
A 

M n'-l 



Pa, q, V, vb] k[ + ,---, k'J] K±, K'±; ^) 



/iVl " — ^ 
n dk[ - j(">--<' (p^, g, r^, k[,..., k'J . (46) 
i=i 

We now want to set fc^ = in order to put Eq. (jl^ into the form of Eq. (001) • To that end, we 
employ an identity for ^ (we again suppress the dependence on the color indices for brevity) 

JA\PA,q,V,VB;k[,...,k'^,) = J^"^ (pA,q,V,VB;k[-^ = 0, = 0,k[-,...,k'^7,k[^,...,k'^,^ 



n'-l 



k - 



i=l Oh ^ 



Substituting the first term of Eq. (|T7j) into Eq. we recognize the definition for Va, Eq. PT|) . 

We have shown in Sec. 14.31 that the contributions from the terms proportional to dJ^^/dlf in 
Eq. (jTrf) can be expressed as soft-loop integrals of some \ again with fewer loops than in 
\ When we substitute this into Eq. (j46p we may express the resulting contribution in terms 
of integrals which have the form of Eq. ()42|). We can now repeat all the steps mentioned so far, 
with this new integral. By this iterative procedure we can transfer the k~ integrals in Eq. (jl^ to 



J^"^ and also set k^ = inside J^"^ In a similar manner, we can analyze the integrals in Eq. 
(|27p. and express them in terms of defined in Eq. (j41j) . This algorithm, indeed, leads from 
the first factorized form of the considered amphtude, Eq. (j27j), to the second factorized form, Eq. 
iOl). 



5.2 Evolution equation 

We have now collected all the ingredients necessary to derive the evolution equations for quantities 
defined in Eq. pij) . Consider F^\ We aim to find an expression for p'^dT^^'' / dp\. As discussed in 
Sec. 14.21 this will enable us to resum the large logarithms of In(pJ^) and eventually the logarithms 
of ln(s/|t|). According to Eq. (HH), in order to find p^dT^^^ / dp\, we need to study p\dJ^J^'' / dp\. 
Using the identities Pa- vb = '^^aCb, Pa ■ h = 2^A^i, where = k~ri~^ and C,a, (b are defined in 
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Eq. (jSH), we conclude that 



(4J 



From this structure, using the chain rule, we derive the following relation satisfied by which 
generalizes Eq. to with arbitrary number of external gluons. 



Pa 



dJ 



(n) 



dp], 



(n) 



n-1 



dJ 



(n) 



i=l 



dk: 



+ c 



dJ 



(n) 



B- 



dC 



B 



(49) 



Now, we integrate both sides of Eq. 



over 



n;:; (x 



the definition for , Eq. (jUj), the left hand side is nothing else but p\dT^2 / ^p\. The first term 



and set all k^ = 0. Then, using 

in) 



on the right hand side of Eq. (pUj) is simply —ri°'dT^2^/dri°'. Noting that (sdJ]^^ /dC,B = n J 



An) 



An) 



the last term gives simply nF^^ For the middle term, we use integration by parts 



n-l / nM \ n-1 o j{n) n-1 / „M s ra-1 

n(/.,*7)E*r|^ = n(/„*7)E 

VJ_M / Cfci . 1 \J-M J 



dk; 



A ) 



J 



in) 



n-l f-M /n-l 



i=l 



M 



M 



4Hkr = +M,...) + jr{k- = -M,... 



in). 



(n-l)F?).(50) 



Combining the partial results, Eqs. and (fSUj). we obtain the following evolution equation 



Pa 



n-l pM (n-l 



1=1 



M 



M 



j'2\k- = +M,...) + J2>{k; = -M,...) 



in), 



_1_ r^n) 



in) 



(51) 



The jet function J^""* in the first term of Eq. (j3T|) is evaluated at {kf = 0}f^i and the k~s are 
integrated over for j = 1, . . . , n — 1 and j ^ i. The first term in Eq. ()51|) can be analyzed using 
the K-G decomposition for gluon i since the k~ is evaluated at the scale M ~ \/[t\. The outcome 
of the last term in Eq. (jKT|) has been determined in Sec. 14. 2| Fig. |3]^. As a result we have all the 
tools necessary to determine the asymptotic behavior of the high energy amplitude for process 
(0). To demonstrate this, we will rewrite Eq. ()5ip into the form where on the right hand side 
there will be a sum of terms involving F^ ■'s convoluted with functions which do not depend on 
p\. Let us proceed term by term. 

Again, the K-G decomposition applies to the first term in Eq. ()5H) because the external 
momenta are fixed with k~ = ±M. Using the factorization of a gluon given in Eq. ()3U|) it 



^Strictly speaking we have analyzed fi" d j')^^ / d 1]°' , but because of the relationship between 



(n) 



and r^'' given 



by Eq. 141|l . once we know f)'^ d j'^^ / d ri°' we also know fj^dTj^^/dri" 



(n) ^ 
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is clear that the contributions from the K gluons cancel for jji^^s evaluated at k~ = +M and 
k~ = —M. Hence only the G gluon contribution survives in this term. Its most general form is 
shown in Fig. ^jp. Before writing it down let us introduce the following notation. For a set of 
indices {1, 2, ... , n}\{i} consider all the possible subsets of this set, with 1, 2, ... , (n—l) number of 
elements. Let us denote a given subset by n, its complementary subset vf, the number of elements 
in this subset as tZtt and in its complementary as = (n — 1) — tZtt. With this notation, we can 
write the ith contribution to the first term in Eq. ()5H) in the form 



„ N-l 

jin)a,...a„ (A;r = +M, . . .) + jj^ ■ (k^ = - M ,...)= J n 



\ (2vr) 



D 



^a^a^^.'^ai^^ 61... fe]v (^i ~ +M, fcj^ , . . . , k^^^ ; k^' — 0, k^^ — 0, . . . , kj'^^ — 0; ki±, ki-^±, . . . ,k 



{nff+N)aj ...aj bi...bisr / _ _ . . 

h,...jN;PA,q,v) + iK -M). (52) 

In Eq. (1221), the summation over repeated indices is understood. We sum over all possible subsets 
TT. In other words, we sum over all possible attachments of external gluons to jet function Ja and 
to the soft function 5*. The elements of a given set vr are denoted ii,i2, . . . , in„. The elements of a 
complementary set vr are labeled zi, ^2, ■ • ■ , ^n*- The number of gluons connecting 5* and J^""""*"^^ 
is N. 

Following the procedure described in Sec. 15.11 with Ra in Eq. replaced by S in Eq. (jS21), 
we can express the contribution from a G gluon in the first term of Eq. (|HT|l in a form 

1 „j\,f /n— 1 



Y.J \^^^'^] ^ [jf''""''^(K = +M,...) + jf'''-''-{kr = -M,...] 

/m 
n d^-'hi. • • • , ^n±, . . . , Im^-^ g, r/; M) 

■III 'j=l 

X r^^^'^-''"(p^,g,r7;/i^,...,W;M). (53) 

The function /C("'™') does not contain any dependence on pA- It can contain delta functions setting 
some of the color indices 6j, as well as transverse momenta of equal to color indices and 
transverse momenta fcj_L of F^\ 

Next we turn our attention to the last term appearing in Eq. ()51|) . The contribution to this 
term has been depicted graphically in Fig. |3] Consider the term in Fig. It can be written in 
a form 

/M /n-l \ f-M /"-I \ 

/c^ = 0, A;^ = 0, . . . , kf^^ = 0; ki ±, ki^ ±, . . . , ki^^ _|_; Z = fcj + fcj^ + . . . + ki^^ 

X ^ (fcjj , • • • , kj^_ ; = 0, . . . , = 0; /Ci^ _L, . . . , /Ci„_ _l; 

l = ki + ki^ + ... + ki ■,pA,q,v)- (54) 
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In Eq. (|54p . we have used the same notation as in Eq. (|52j) . Momentum / connects S' with 
jjn^+i)_ Following the same procedure as in Sec. 15.11 with Ra appearing in Eq. (j42p replaced by 
S' introduced in Eq. (j54j) . we can express this contribution in a form given by Eq. (j53p with a 
different kernel /C^"'"^). 

The contribution from Fig. HJd can be written 

Ja = 0, = 0;/4^_L, ...,/Ci„_^;fc,/^ (55) 

The flow of momenta k and / is exhibited in Fig. ^p. The momentum k flows through the boxed 
vertex and the ghost line shown in Fig. EJd which forces this momentum to be soft, so that lines 
k and / are part of the function S" . Since the line with momentum k is soft, then all gluons 
attaching to in Eq. ^ are soft and we can again apply the procedure described in Sec. 

l5.1l to bring the contribution in Fig. 0]d into the form given by Eq. (j^Hj) with a different kernel, of 
course. 

In summary, we have demonstrated that all the terms on the right hand side of Eq. (j^Tj) can be 
put into the form given by Eq. This indicates that Eq. (jHH), indeed, describes the evolution 

of F^'' in hip\ since it can be written as 

(^^ 9p+ - 1^ ry''"''"(pA,g,^/;A;i±,...,A;„x) = 

m j=l 

xT^^^''^-'-{pA,q,vJi±,...JmL). (56) 

The kernels /C("''") do not depend on p\. As indicated above, they can contain delta functions 
setting some of the color indices bi, as well as transverse momenta li± of F^'' equal to color indices 
and transverse momenta ki± of F^\ The systematic use of this evolution equation enables us 
to resum large logarithms \ii{p\) at arbitrary level of logarithmic accuracy. Analogous equation 
is satisfied by F^. It resums logarithms of In(p^). 



5.3 Counting the number of logarithms 

Having derived the evolution equations for F^'', Eqs. (jSH) and it does not take too much 
effort to show that at r-loop order the amplitude contains at most r powers of ln(s/|t|). We follow 
the method of Ref. jlj. We have argued in Sec. 15.11 that the power of ln(s/|t|) in the overall 
amplitude corresponds to the power of In(pJ^) in F^^ So we have to demonstrate that at r-loop 
order F^''^\ where F^'^'' represents a contribution to F^^ at r-loop level, does not contain more 
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than r logarithms of \n{p\). We prove this statement by induction. First of all, the tree level 
contribution to r^'°^ is proportional to the expression 

/M /n-l \ n-l / ^ \ 

/ I., .J ,=1 - S.1 fc.)' + 'I \,.„ ) 

where t'^'^ s are the generators of the S'f/(3) algebra in the fundamental representation. The sum 
over indicates that we sum over all possible insertions of the external soft gluons. 

Eq. (|Ffj) is evaluated at {kf = 0}[L^. Expanding the denominators in Eq. ^7\i we obtain the 
expression —2p~^{k^^ + . . . + k^^) — {ki^ + . . . + ki.)\ + ie. We see that the poles in k^ planes are 
not pinched and therefore the k~ integrals cannot produce In(pj^) enhancements. 

Next we assume that the statement is true at r-loop order, and show that it then also holds 
at (r + l)-loop level. To this end we consider the evolution equation, Eq. ()51|). and examine 
{p\ d/dp\ — l) r^''^^^''. Its contribution is given by the first and the third term on the right hand 
side of Eq. fl51|) . As already mentioned, the first term in Eq. (jSip can be analyzed using K-G 
decomposition. The contributions from the K terms cancel each other while the contribution from 
the G gluons are given by the kind of diagram shown in Fig. The latter, however, can be 
written as a sum of soft loop integrals over ^ with r' < r, since we loose at least one loop 
in the original J^'^^^^ due to the soft momentum integration. This is demonstrated in Eq. (j52p . 
Following the procedure described in Sec. 15.11 we may express these contributions as transverse 
momentum integrals of some F^''^'', see Eq. (j53j) . These contain at most r' < r logarithms of 
In(pJ^). The contribution from the third term in the evolution equation, Eq. ()51|). is given by the 
diagrams depicted in Fig. EJ These are again soft loop integrals of some J^" ^ with r' < r, and 
they can be expressed as transverse momentum integrals of F^ \ see Eqs. ()54p and which 
have, therefore, at most r logarithms of In(p^). Since both terms on the right hand side of Eq. 
dnH) have at most r logarithms of In(p^), then also p\dT'^2''^^^'^ / d p\ has at most r logarithms of 
ln(pj[) at (r + l)-loop level. This immediately shows that p^'^^^'' itself cannot have more than 
(r + 1) logarithms of In(p^) at (r + l)-loop level. 

This result enables us to formally classify the types of diagrams which contribute to the 
amplitude at the fc-th nonleading logarithm level. As has been shown in Sec. 15. ![ we can write an 
arbitrary contribution to the amplitude for process (P) in the Regge limit in the second factorized 
form given by Eq. (pOj) . Consider an r-loop contribution to the amplitude and let L^, Lb and Ls 
be the number of loops contained in F^, F^ and S. Since F^ (F^) can contain La {Lb) number 
of logarithms of p\ (j9^) at most, the maximum number of logarithms, A^'maxLog, "we can get is 

iVmaxLog = T - Lg. (58) 

This indicates that when evaluating the amplitude at the fc-th nonleading approximation, we need 
to consider diagrams where 1, 2, . . . , (fc + 1) soft gluons are exchanged between the jet functions 
J A and Jb- 
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5.4 Solution of the evolution equations 

Having obtained the evolution equations, Eqs. and we discuss how to construct their 
solution. Our starting point is Eq. (j^Uj) . In shorthand notation it reads 



P\ ri"^^ = E E ^^""'""^'^ ® ^a''\ (59) 
Pa ^/^o n' 

at r-loop level. Indices n and ra', besides denoting the number of external gluons of the jet function, 
also label the transverse momenta and the color indices of these gluons. The symbol ® in Eq. 
fl59|l denotes convolution over the transverse momenta and the color indices. Note that Eq. 
holds for Va with the overall factor p\ divided out {Va = '^a/Va)- have proved, in Sec. 15.31 
that r^'^"* can contain at most r logarithms of In(pJ^) at r-loop level. Therefore the most general 
expansion for T a is 



r 



j=0 

If we want to know F^''^'' at N'^LL accuracy (fc = is LL, = 1 is NLL, etc.), we need to find all 
c'p^^ such that r ~ j < k. The coefficients c^''^'^ in Eq. ()6U|) depend on the transverse momenta 

and the color indices of the external gluons. Using the expansion for r^"''^) and T^^' '"'^ Eq. dnUD, 
in Eq. (|59p and comparing the coefficients with the same power of In(p^), we obtain the recursive 
relation satisfied by the coefficients c'p^^ 

r—l n+r—r' 

^^in,r)^ E E (61) 
r'=j — l n'=l 

In Eq. fl6H) . we have used that, in general, l<n'<n + r — r'. 

We now show that Eq. (jUT|) enables us to determine all the relevant coefficients c'p"'^ of 
order by order in perturbation theory at arbitrary logarithmic accuracy. We start at LL, k = 0, 
and consider n = 1. At r-loop level we need to find the coefficient ci^'^\ It can be expressed in 
terms of lower loop coefficients using Eq. (jblj) and setting j = r and n = 1 

r4^'^) = ^^(i.<i)®ci!:r'^. (62) 

n'=l 



In Sec. I6.1l we will prove that the one loop kernel satisfies /C^^'^'^^ = 0, Eq. ()72|) . This implies that 
in Eq. ()62j) the coefficient Cr^'^'* is expressed in terms of lower loop coefficient c^^-T^^ and hence, 
we can construct the coefficients at arbitrary loop level once we compute Cq^'^"*, the coefficient 
corresponding to the tree level jet function F^'"''. 

Next we construct all F^"* for ri > 1 at LL accuracy. Let us assume that we know all Cr" '^^ for 
all r and for n' < n. We apply Eq. ()61|) for j = r 



r Cj. 



n+l 



n'=l 
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In Sec. 16.21 we will show that the evolution kernel in Eq. (j63|) obeys /C^"'"''^) = 6(n — n') 

Eq. (jlUSp . where 6{n — n') is the step function. This implies that the sum over n' in Eq. (j63j) 

terminates at n' = n. Isolating this term in Eq. (jUHj) . we can write 

n-l 
n'=l 

So after we calculate the tree level coefficient Cq"'"'' , we can construct all the coefficients Cr"''^'' using 
Eq. (j64p order by order in perturbation theory, since according to the assumption we know cf" '^^ 
for all r and for n' < n. This proves that we can construct the jet functions at LL, /c = 0, for all 
n to all loops. 

We now assume that we have constructed all the jet functions at the N'^LL accuracy for a 
given > and we will show that we can determine all the jet functions at the N'^+^LL level. We 
start with n = 1. Using Eq. with n = 1, j = r — (/c + 1), isolating the term with r' = r — 1 
in the sum over r' and using A^(i'"''i) = S^^^, }C^^'^'^\ we arrive at 

ir-k-i) et, = }c^'''-^'^ ® eri + £ 'y: ® (65) 

r'=r—k—2 n'=l 

After we evaluate the coefficient c^''^'^^^ (impact factor), Eq. implies that we can calculate 
the coefficients c'r-k-i order by order in perturbation theory, because, according to the induction 
assumption, we know all the coefficients c^" jr_2 since they are at most N'^LL. Once the coefficients 
of r^'' are determined at N^"'"^LL level, we assume that we know all the coefficients of ^s for 
n' < n. We want to show that we can now construct all the coefficients for F^^^ at N^'^^LL 
accuracy. First we need to calculate c^^'''~^^\ Then we use Eq. to express the coefficient 

c"J!fk-i, isolating the terms with r' = r — 1 and n' = n, as 

n— 1 
n'=l 

1 — 2 n+r—r' 

+ y: e (66) 

r'=r—k—2 n'=l 

The terms appearing in the sum over r' in Eq. ()66|) are known according to the assumptions 
since for them r' — {r — k — 2) < k. We also know, according to the induction assumptions, the 
contributions to the second term of Eq. ()66|) . since they have n' < n. Therefore, we can construct 
c^^-k-i order by order in perturbation theory. This finishes our proof that we can determine the 
high energy behavior of F^^ at arbitrary logarithmic accuracy. Note that to any fixed accuracy 
only a finite number of fixed-order calculations of kernels and coefficients Cq"'*^^ must be carried 
out. In a similar way we can construct a solution for F^\ 

Once we know the high energy behavior for F^^ and F^\ then the second factorized form, 
Eq. pni) . implies that we also know the high energy behavior for the overall amplitude. Because 
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a) b) 

Figure 6: Diagrams contributing to the amplitude at NLL approximation: factorized one gluon 
exchange diagram (a) and nonfactorized two gluon exchange diagram (b). 

a jet function F^"-* is always associated with at least n — 1 soft loop momentum integrals in the 
amplitude, we infer from Eq. ()58j) that if we want to know this amplitude at N^LL accuracy, it 
is sufficient to know F^^"^ (F^^^) at N-^+^-^LL (N^+^-'^LL) level for n < K + 1 {m < K + 1). 
We note, however, that to construct these functions according to the algorithm above, it may be 
necessary to go to slightly larger, although always finite, values of n and m. Let us describe how 
this comes about, starting with the basic recursion relations for coefficients, Eq. ()6ip . 

We assume that for fixed n on the left-hand side of Eq. (j6T|l . the logarithmic accuracy k is 
bounded by the value necessary to determine the overall amplitude to i^th nonleading logarithm: 
k = r— j<K + l— n, which we may rewrite as n + r — {K + 1) < j < r. On the right-hand side 
of Eq. ()61|) we encounter the coefficients of the jet functions with n' external lines, satisfying the 
inequality n'<n + r — r'<n + r — {j — 1). Combining these two inequalities, we immediately 
obtain that n' < K + 2. Then, for any given number of external gluons n' on the right-hand side, 
we encounter a level of logarithmic accuracy k' = r' — {j — 1) < n + r — n' — {j — 1) < K -\-2 — n'. 
This reasoning indicates that, in general, we will need all fJ^"'^ (fJ"'^) at N^+2-"'LL (N^+^-^'LL) 
level for n' < K + 2 (m' < K -\- 2), when evaluating the amplitude at N^LL accuracy. We note 
that for fermion exchange in QED it was shown in Ref. [4J that only contributions with n' < K-\-l 
are nonzero, but for QCD, two-loop calculations appear to indicate, Ref. that QCD requires 
the full range of n' identified above, starting at NLL. 

6 High energy behavior of the amplitude 

In the previous sections we have developed the general formalism for obtaining the high-energy 
behavior of the scattering amplitude for process (P) at arbitrary logarithmic accuracy. In the 
following subsections we apply these techniques to study this amplitude at LL and NLL level. 

6.1 Amplitude at LL 

According to Eq. ()58p. the amplitude at LL comes solely from the factorized diagram shown in 
Fig. inti, but without any gluon self-energy corrections. The jet A, containing lines moving in 
the plus direction, and jet B, consisting of lines moving in the minus direction, interact via the 
exchange of a single soft gluon. This gluon couples to jet A via the — component of its polarization 
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Figure 7: Diagrammatic representation of the evolution equation for jet J^^"* at LL. 

and to jet B via the + component of its polarization. Since v^Nap{q,rf) = 1, we can write at 
LL ^ 

^8^8 = —^JA^''{PA,q,v) JB^''{PB,q,v), (67) 

where feg is the color basis vector corresponding to the octet exchange, defined in Eq. (jS)). Using 
s = 2p^p^, the logarithmic derivative of the amplitude can be expressed as 

dAs^ __ia/^ 1 (i)„a7^ 

In Sec. 14. 2[ Eq. (jHHj) . we have derived an evolution equation resumming In(p^) in J^^\ We note 
that = T^}\ and that ^ is a special case of the evolution equation ()5ip . The diagrammatic 
representation of the first term on the far right hand side of Eq. (jHHj) . which follows from Fig. 
m in the case when we have one external soft gluon attached to a jet function, is given by the 
diagrams in Fig. [7| Diagram in Fig. [7^ corresponds to Fig. EJd and the diagrams in Figs. [Tjp and 
[7\: correspond to Fig. for n = 1. The diagrams in Figs. [Tjp and [3; are in the factorized form, 
while the one in Fig. [7^ is not. 

As discussed in Sec. 14. 2[ power counting shows that the loop momentum k in Fig. [7^ must be 
soft. This implies that we can make the following approximations. First, since at LL all internal 
lines of the jet A are collinear to the + direction, we can neglect the dependence of i.e. 
we may set = inside J\ . Also, we can pick the plus components of the vertices where the 

(2) 

soft gluons attach to the jet J\ . A short calculation, which uses the Feynman rules for special 
lines and vertices listed in Appendix El gives the contribution to Fig. in a form 

X vl jf})l^ {pa, q, V, = 0, k' , k^) , (69) 

where we have defined Qs = gs^^'^■ Using Eqs. (fTT|l and for the components of the gluon 
propagator and the boxed vertex, respectively, it is easy to see that in the Coulomb (Glauber) 
region, k~ <^ k^ k±, the integrand in Eq. (jU^ becomes an antisymmetric function of k~^ and 
that therefore the integration over k~^ vanishes in this region. 
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a) b) c) 

Figure 8: Diagrams determining the contributions to the gluon trajectory at the order 



In the soft region, where all the components of soft momenta are of the same size v — t, we 
can use the K-G decomposition for the soft gluon with momentum k attached to J^^"*. At LL, 
however, there cannot be any soft internal lines in j'^^^ in Eq. since, as discussed in Sec. 15.31 
only integrals over collinear momenta can produce powers of InpJ^. Therefore, at LL, only the 
K gluon contributes, because the G gluon must be attached to a soft line. The K gluon can be 
decoupled from the rest of the jet jjf"* using the Ward identities, Eq. (j30|) . Their application in 
Eq. (jnni) gives 

X J^'-lp^,?,!)). (70) 

We have used the identity fad, fdcb = 5ad = Ca Sad in Eq. (fTOj) . Eq. (fTUj) now gives a factorized 
form for Fig. [7^. Since the contributions in Figs. [Tjp and [Tt: are already in the factorized form, 
we can immediately infer that the gluon reggeizes at LL. Combining the terms from Fig. [7|in Eq. 
dnSl), we obtain the evolution equation at leading logarithm 



d 



PA-g-^JA^''iPA,q,v) = a{t) j'^^^{pA,q,v)- (71) 

Using the notation for evolution kernels introduced in Sec. 15. 4[ Eq. (|7ip implies that 

/C(i'2;i) = o_ (72) 

In Eq. dZH) 

a{t) = 1 + + + (73) 

is the gluon trajectory up to the order a^, and a^a\t), al^\t) and a''c\t) are its contributions 
given in Figs. ^ - |Hfc, respectively, 

d^k 1_ 

(27r)^ k^{k-qyk -kvA-k 

a^^\t) ^ i-glCA j (0? fc2(fc_g)2fc.fc ^"(^)^""(^) <N;{q-k)V,Uk,-q,q-k)vl 
^\^) - -^-^ICa j ,.,.,(^_\).(^-_,) <K,m) i^s ■ k) . (74) 



a«(t) = -^-glCAt j—^-— —-^ -v'^N,,ik)S'^ik)v^Mq-k)v''A, 
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In Eq. (fTij). V^i3^{k, —q, q — k) stands for the momentum part of the three-point gluon vertex. 
After contracting the tensor structures in Eq. (f7i|) . using the exphcit form for V^i3^, Va, Vb, S'^ 
(Eq. (jMjl ) and for the components of the gluon propagator, Eq. (fTT|. we obtain for a^^bci^)^ 



t f d^k [klko + k^ks][{k-qy -{k-qf_ 



2 J (27r)^ {ko + ks) P {k - g)2 (A; ■ k^ {k - q) ■ [k - q) 

1 r d^k 1 



[fcl (A; - q^ + 2k'^ kj {k - g)^ ■ + 2kl k] k± ■ {k - q)± + 2kl k\ {k - qf ' 



n 5 



Next, we perform the k^ and A;^ integrals in Eq. (|75p. For ai^''(t), these integrals are UV/IR finite. 
However in the case of cx^^l{t), the k^ integral is linearly UV divergent. In order to regularize 
this energy integral, we invoke split dimensional regularization introduced in Ref. [HZj- The 
idea is to regularize separately the energy and the spatial momentum integrals, i.e. to write 
d^/cg d^^k^d^'^k for Euclidean loop momenta kE- The dimensions Di and D2 are given by 
Di = 1 — 2ei and D2 = 3 — 2e2, with ej —>■ 0+ for j = 1, 2. Since the energy integral for a^c\t) 
is scaleless, it vanishes in this split dimensional regularization. The energy integrals in al^l{t) are 
st raight forward . 

All the k^ integrals can be expressed as derivatives with respect to k'j_ and/or {k — q)\ of a 
single integral 



Jo ^Jkl + a2 {kl + 62) hy/h^ ~ a? \ a J 

The result of these integrations over k"^ is 

riw N 2e^ fd^-^k^f^,,, , „ kl 2(k-q)l-3kl 

a«(i) = a^f^'^CAt —-r^ili\k^\,\k±-q^\)j-, ,21. + — — 



(27r)^-2 V " "' ^"'^ [{k-q)l-kl]^ kllik - q)l - kl]^ 



a,{t) = -^sf^'CAt J j^,^I{\k,\,\k,-q,\)^^^_^^,^_^,^^^ [ik-q)l-kl]\ 

a^}\t) = 0. (77) 

Combining the results of Eq. (fTTjl and Eq. (|7H|l . we obtain the standard expression for the gluon 
trajectory at LL 

a{t) = 1 + C^a,^- I "^^^^l^ . (78) 

We can now simply solve the evolution equation (jUH|) . to derive the factorized (reggeized) form for 
the amplitude in the color octet 

As{s,t,as) = s''^'^As{t,as). (79) 

The amplitude factorizes into the universal factor s"^^\ which is common for all processes involving 
two partons in the initial and final state and dominated by the gluon exchange, and the part A^, 
the so-called impact factor, which is specific to the process under consideration. 
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Figure 9: Expansion of the one gluon exchange amphtude at NLL using the tuhp garden formal- 
ism. 

6.2 Amplitude at NLL 

At NLL level the contribution to the amplitude comes from both the one gluon exchange diagram, 
Fig. inti, and from the two gluon exchange diagram, Fig. IHb- At this level, both singlet and octet 
color exchange are possible in the latter. Including the self-energy corrections to the propagator 
of the exchanged gluon (taking into account the corresponding counter-terms), we can write the 
contribution from the diagram in Fig. IHt as follows, 

A^'^ = -] J['J)1 (pa, g, V) (^N-^q, v) + ]v% < U,,{q, ^) t;^ v^^ (p^, g, r^), (80) 

where Yi^y{q,rj) stands for the one loop gluon self-energy. We now put this contribution into the 
first factorized form, Eq. ()27|). isolating the plus polarization for jet A, and the minus polarization 
for jet B. At NLL in the amplitude, we need the soft function S'^^'^\ Eq. (j27j) with n = m = 1, to 
accuracy 0{as). Using the tulip-garden formalism described in Appendix [01 the contribution to 
the first term on the right hand side of Eq. (jHn|) is given by the subtractions shown in Fig. El In 
accordance with the notation adopted in Appendix O the dashed lines indicate that we have made 
soft approximations on gluons that are cut by them. A dashed line cutting a gluon attached to jet 
A[B) means that the gluon is attached to the corresponding jet through minus(plus) component 
of its polarization. Since = in the Regge limit, Eq. (0), we have A^^^(g) = g^^. This implies 
that the contributions between the diagrams in Fig. and in Fig. Eli as well as between the 
diagrams in Fig. Et and in Fig. EF cancel each other. Therefore only the zeroth-order soft function 
diagram in Fig. Eb survives in the factorized form, Eq. ()27p. 

For the two gluon exchange. Fig. (Hb; we only need the lowest order soft function at NLL in 
the amplitude (and LL in singlet exchange). The expression for the two gluon exchange diagram 
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in Fig. takes the form, Eq. (|77j). 



where ^(/c) is given by 

We have suppressed the dependence of the functions appearing in Eq. (|HT|) on other arguments 
for brevity. At NLL accuracy we are entitled to pick the plus Lorentz indices for jet function 
and the minus indices for jet function Jb only. We can also set /c"*" = in and = in Jb 
since all the loop momenta inside the jets are collinear. Eq. 1)811) represents the first factorized 
form, Eq. (j2Zj), for the amplitude A^'^\ 

Next, we follow the procedure described in Sec. 15.11 to bring the amplitude into the second 
factorized form, Eq. (jlO)). We employ an identity based on Eq. f|44|l . for the function S{k) defined 
in Eq. ^ 

s{k+,k-) = s{k+ = o,k- = o)e{M -\k+\)e{M -\k-\) 

+ [S{k+, fc- = 0) - s{k+ = o,k' = 0) e{M - \k+\)] e{M - \k~\) 

+ [S{k+ = 0, k~) - S{k+ = 0,k- = 0) e{M - |A;~|)] e{M - \k+\) 

+ [{S{k+, k-) - s{k+, k- = 0) e{M - \k-\)\ - 

{S{k^ = 0, k-) - S{k+ = 0,k- = 0) e{M - \k-\)} 0{M - \k+\)]. 

(83) 

The contribution from the first term in Eq. (|83|) gives immediately the second factorized form 



m 



with T]^ and defined in Eq. ()4ip for n 

We now discuss the rest of the terms in Eq. (jHSl), which can be analyzed using the K-G 
decomposition, since, by construction, there is no contribution from the Glauber region. At the 
current accuracy only the ii'-gluon contributes. After substituting the second term of Eq. ()83|) into 
Eq. (jHH), we can factor the gluon with momentum k from jet However, it is easy to verify, 
using the definitions for K and G gluons, Eq. ()29|) . the Ward identities, Eq. ()3Up . and the explicit 
components of the gluon propagator, Eq. (jlip . that the k~^ integral is over an antisymmetric 
function. As a result, this contribution vanishes. In a similar fashion, the contribution from the 
third term in Eq. ()83p. after used in Eq. ()ST|) . vanishes, since now we can factor the soft gluon 

(2) 

with momentum k from jet J\ and the k~ integral is over an antisymmetric function. 

In the case of the last term in Eq. ()83|) . after used in Eq. ()HT|) . we can factor the soft gluon 
with momentum k from both jets J\ and . The integrals of the soft function S{k) over fc"*" 
and k~ are then 

^2 f^^ dk+ dk- 



S{k^, q- M) ^ Ca ^ j ^^-^ S{k^, k^, q). (84) 

As usually, we leave the transverse momentum integral undone. The l/fc"*" and 1/A;^ in the 
integral above are given by the Principal Value prescription because there is no contribution from 
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the Glauber region. Since the amphtude is independent on the choice of scale M, we can evaluate 
it at arbitrary scale. We choose to work in the limit M ^ 0. In this limit the contribution to 
the integral comes from the imaginary parts of the gluon propagators in Eq. (jH^ . —iT[5{k'^) and 
—m5{{k — qY). The integration is then trivial and Eq. ()84|) becomes 

. lim S(*.,,;M) . -C, f j^^J-^. (85) 

Combining the partial results of the analysis described above in Eq. (jHT|) . we arrive at the second 
factorized form for the double gluon exchange amplitude, Fig. IHb, 

A''' = I "^(^x) ^ sik-^ = o,k~ = 0, k^) rg) ^\k^) 

^-^^ r« "(p^, q) Sik^, q) r« \pA, q). (86) 

Using Eq. (|Hn|l for A^^'^ and Eq. (|HF)jl for A^^^ , we obtain the amplitude for the process ((H) at NLL 
accuracy 

= -irw-ip^,,,,,) ('i + in+ _(,,„) + 1 £,<"(«)) Tf{ps,q„)) 

- /g^^r-',*jgl,,^r-'(*.). (87) 

In Eq. (jHTjl . we have used the explicit form for ^(fc^ = 0, = 0,fc_L), which can be easily 
identified from Eq. (jH^ . We have also used the integral representation of the gluon trajectory 
given in Eq. (fTSj) . 

In order to determine the high energy behavior of the amplitude in Eq. (I87p . we need to 
examine the high energy behavior of or at NLL and the evolution of F^'* or T^"^ at LL. 
In this paper, we restrict the discussion of evolution equations to LL level, and hence we analyze 
the behavior of V\ only. We will address the study of NLL jet evolution, and gluon reggeization 



at this level, elsewhere 

We use the evolution equation given by Eq. ((^T| in order to determine the LL dependence of 
F^-* on In(pJ^). In our special case of the two gluon exchange amplitude, it reads 
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p+^-l]vf''' = M J'^>"\k- = +M,k+ = Q,k^) + J'^"'\k- = -M,k+ = Q,k, 



7'(2)a^'/L-_ I /\/r U+ - Pi U \ I 7(2)a6/L- 



d p(2) ab 



The first term in Eq. ()88j) can be analyzed using the K-G decomposition. The contributions 
from the /T-gluon cancel between the J\ {k~ = +M) and J\ {k~ = —M). The contributions 
from the G gluon, which we now discuss, are shown in Figs. ITIlk andllOb. 

Since the gluon with momentum g — in Fig. ITUk cannot be in the Glauber region, we can use 
K-G decomposition on it. The K part factors from while the G part does not contribute at 
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LL. After factoring out the gluon with momentum q — k and performing the approximations on 
the jet function the contribution to Fig. ITUk for k~ = +M is 

Fig.m = -^9lfaecfdel>l^ J 5i = 0, fc" = +M, fc^, /) ^''(/+ = 0, /" , Z^) , (89) 

where we have defined 

Next we follow the established procedure. First, we write 

s^{k,i) = Si{k,i- = o)e{M-\i-\) + [Si{k,i) - Si{k,i- = o)e{M-\r\)] . (91) 

When we use the second term of Eq. (|9ip in Eq. (j89p . we can factor the gluon with momentum 

(2) 

/ from J)^ . Since the resulting integrand is an antisymmetric function under the simultaneous 
transformation M —M, the contributions on the right hand side of Eq. (jHHj) 

evaluated for k~ = +M and k~ = —M cancel each other. Therefore we can write, using Eq. (PT|) 
in Eq. dHl, 

Fig.Unt = -tg^faecfdebj^ J dl+ Si{k+ = , k' = +M , k^, = , h) T^^^ (h) 

+ (92) 

where by dots we mean the term which is canceled after we take into account the contributions 
to both J^^\k' = +M) and J^\k' = —M) on the right hand side of Eq. 
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Next, we perform the integral in Eq. (j92j) . As we have already mentioned above, since the 
final result does not depend on the scale M, we can choose arbitrary value of M. We have chosen 
to perform the calculation in the limit M —>■ 0. Then the only nonvanishing contribution comes 
from the imaginary part of the propagator !/[(/ — + ie], —iir 6{2Ml~^ + ~ this term 

the integration is trivial and we obtain 

M (Fig. m) = faeche, j ^t^^ T^^ {I ^) + . . . , (93) 

which gives an M-independent contribution to the right hand side of Eq. (jEHI)- 

We follow the same steps when dealing with the diagram in Fig. ITUb . whose soft subdiagram 
is given by 

S2{k,l) = — ^ -, -^Vf,p^{l,-k,k - I) (g^^ -Tj T^V^s-[q- 1,1 - k,k - q). 



(94) 

First we use the identity for The contribution due to the second term in Eq. vanishes, 

(2) 

after the gluon with momentum / has been factored from J\ , due to the antisymmetry of the 
integrand. Hence again, as in the case discussed above, only the term given by 5*2 (/~ = 0,l~^,lj_,k) 
contributes. In the limit M — >■ 0, the contribution comes from the imaginary part of the same 
denominator as in the case of Fig. ITUk . The result is 



M(Fig.[inb)= - asfaecfdeb J 



d^-2/ 



X 



(27r)^-2 II (/ - q)l {k - 1)1 
{kill - k^ ■ IJl - k^ ■ qJl - klh ■ q± + 2k^ ■ IJ^ ■ gx) x vf^^lj 



+ .... (95) 
Combining the results of Eqs. (j93|l and (j95|l . we obtain the expression for the surface term in Eq. 



M 



jf ^\k- = +M, k+ = 0, ki_) + jf "'{k- = -M, k+ = 0, k, 



faecfbed 



J {2n)^-^[ll{k-l)l + {l-q)l{k-l)l ll{q-l)l) ""'^ ^^'^ 

Next, we analyze the contributions to the term fj^d/dri^T^^^ in the evolution equation ()88|) . 
The contributing diagrams are shown in Figs. ITUb - [TUF . Note that for every diagram in Figs. ITUb 
- [TUF . we have also diagrams when a loop containing the boxed vertex is attached to the external 
gluon with momentum k, instead of to the external gluon with momentum q — k. 

In Fig. ITUb . we have to consider all the possible insertions of external gluons with momenta 
k and q — k. We have six possibilities. The contribution shown in Fig. ITUb is proportional to 
(omitting the color factor) 
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a) b) c) 




Figure 11: Contributions to the diagram in Fig. imF when the gluon coming out of the boxed 
vertex is attached to the soft hne (a) and when either or both gluons with momenta k and / are 
K gluons and they are factored from the jet (b - f). 



Since the integrand is an antisymmetric function under k~ — ^ —k~ and — > —l^, the integral in 
Eq. ()97|) vanishes. The same antisymmetry property holds for the remaining five diagrams and 
therefore, there is no contribution from them. 

Let us next focus on the diagram in Fig. llOf . When the gluon with momentum / attaches to 
a soft line inside of the jet jjf"*, the contribution takes the form shown in Fig. ITTk . If it attaches 
to a jet line, its contribution can be written as 

Fig.m = -gsfbcd J -^Ssik+ = 0,k-,k^J) /^^'''\k-^ = 0,k~,k^J+ = OJ-,h), (98) 



with the soft function 



,2 'V^-O SM N-+{q-k-l) 



We use the identity for this soft function S^, obtained from Eq. (jHHj) by the replacement k~^ ^ l~ , 

S3{l-,k-) = Ssil' = o,k- =0)e{M ~\l~\)e{M -\k-\) 

+ [Ssir, k~ = 0)- Ssir = o,k~ = o) e{M - e{M - \k-\) 

+ [S3(/- = 0, k-) - S:i{i- = o,k- = 0) e{M - |A;-|)] e{M - \r\) 
+ [{s^ii-, k-) - s^ir, k- = 0) d{M - \k-\)} 

-{S3{r = o,k-)-S3{i- = o,k- = o)e{M -\k-\)}e{M (loo) 

(3) 

to treat the soft gluons with momenta k and / attached to jet . The contribution from the first 
term in Eq. pOOp . when used in Eq. (j^ . vanishes since the integrand S3{k~^ = k~ = 0, k±, l~ = 
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0,/+,/^) is an antisymmetric function of /+, as can be easily checked using Eqs. (fTTj). and 
(|99p . We can apply the K-G decomposition on the gluon with momentum / when treating the 
second term in Eq. ()100|) used in Eq. (jHEl)- At LL only the K gluon contributes. It can be factored 
from the jet function J\ with the result shown in Figs. ITTb and ITTb . In a similar way we can 
treat the gluon with momentum k in the third term of Eq. ()100|) . After we factor this gluon from 

/n\ ____ ____ 

the jet J\ , we obtain the contributions shown in Figs. ITTH and ITTb. In the case of the last term 
in Eq. pOOj) . we can factor out both soft gluons with momenta k and I from jet J\ . The result 
of this factorization is shown in Fig. lilt '. 

Next, we note that the combination of the diagrams in Figs. ITUti. IlUb and ITTb is the same as 
the result encountered in the analysis of the LL amplitude, Fig. |H| We write 

/M 
dk- (Fig.Cnil + Fig. 1133 + Fig.Eb) = - k^) "''{pa, q, k^). (101) 

-M 

where a^^\q — k) in Eq. ()101|) is given by the diagrams in Fig. |Hlwith an external momentum 
q — k = {0~^, 0~ ,q± — k±). In the case when the gluon coming out of the boxed vertex attaches to 
an external gluon with momentum k, we evaluate the one loop trajectory a^^\k±) in Eq. ()101|) . 

To complete the analysis, we have to discuss the diagrams in Figs. (TTk and ITTb - lilt '. In the 
region ~ /_|_, we can factor the gluon with momentum / from the jet function J^\l^ = 0, 
in the case of the diagram in Fig. ITTb . The resulting k~ and integral is over an antisymmetric 
function of k~ and and therefore it vanishes. So the only contribution comes from the Glauber 
region, where we can set /~ = outside Ja\1^ = 0,Z~,/_l). As above, we perform the and 
k~ integrals in the limit M — > 0. The integrand does not develop a singularity in k~ and/or 
strong enough to compensate for the shrinkage of the integration region J_Ij dk- when M 0. 
Hence the diagram in Fig. ITTb does not contribute in the limit M — 0. In a similar way as for 
the diagram in Fig. ITTb . none of the diagrams in Figs. ITTb - ITTt contribute. The diagrams in 
Figs. ITTb - lllb vanish in the M — ^ limit, while in the case of the diagram in Fig. lilt ' the k~ 
and integral is over an antisymmetric function of k~ and 

At this point we have discussed all the contributions appearing on the right hand side of the 
evolution equation (jHHj) . Combining the partial results given by Eqs. p01|) and in Eq. (jHH|) . 
we arrive at the evolution equation governing the high energy behavior of 

(^P\ - 1^ = iaecfhed J (27^)^-2 ^^^^ '"^ (Pa^ ^' 

, ( kj {k-q)l qj \ 

\ll{k-l)l^ {l-q)l{k-l)l ll{q-l)l) 

+ {a'~'\k^)+a^^\q^-k^))xV^l^''\p\,qM)- (102) 

Projecting out onto the color singlet in Eq. ()102|1 . we immediately recover the celebrated BFKL 
equation [Oj. 
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6.3 Evolution of T^") at LL 



We can now generalize Eq. ()102|) to the case of T]^ . The evolution kernel in this case contains, 
besides a piece diagonal in the number of external gluons, also contributions which relate jet 
functions with different number of external gluons 

p(n) ai...fei...fej...a„/ 4. , , , , x 

n 

1=1 

n-l 

+ (103) 

n'=l 

where ®^ denotes a convolution in transverse momentum space. The last term in Eq. (jl(J3|) 
corresponds to the configurations when one or more external gluons attach to a gluon or a ghost 
lines forming the one loop kernel derived for T)^ . Using the notation of Sec. 15.41 we can write 
Eq. p03|) at r-loop order in a form 

(p\^ - 1) r?'") = J2 ^^"'"''^ ® rj' (104) 

V fA / n'=l 

It corresponds to Eq. (jMj) of Sec. 15.41 when written in terms of the coefficients cl^'^'' introduced in 
Eq. (jUn|) . From Eq. ()104|) we immediately see that the following property of the one loop kernel 

J^{n,n';l) ^ _ ^/^ j^{n,n';l)^ (^05) 

is satisfied. We recall that this step was essential in demonstrating that the set of evolution 
equations, Eq. (jSH), forms a consistent system, refer to the paragraph above Eq. 

The term diagonal in the number of external gluons in Eq. ()103|) coincides with the evolution 
equation derived in Ref. [Sj. Our formalism, besides enabling us to go systematically beyond LL 
accuracy, Ref. indicates that even at LL, in addition to the kernels found in Ref. the 

kernel has contributions which relate jet functions with different number of external gluons. 

7 Conclusions 

We have established a systematic method that shows that it is possible to resum the large loga- 
rithms appearing in the perturbation series of scattering amplitudes for 2^2 partonic processes 
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to arbitrary logarithmic accuracy in the Regge hmit. Up to corrections suppressed by powers of 
|t|/s, the amphtude can be expressed as a sum of convolutions in transverse momentum space 
over soft and jet functions, Eq. (jlHI)- AH the large logarithms are organized in the jet functions, 
Eq. ()41|) . They are resummed using Eqs. (j^T|) and/or (jSHI)- The evolution kernel /C in Eq. (jSUI) is 
a calculable function of its arguments order by order in perturbation theory. This is the central 
result of our analysis. 

As an illustration of the general algorithm we have demonstrated it in an action at NLL for 
the amplitude and LL for the evolution equations. We reserve the study of the NLL evolution, 
which addresses the reggeization of a gluon at NLL, for future work [5^ . 

The derivation of the evolution equations and the procedure for finding the kernels was given 
above in Coulomb gauge. Clearly, it will be useful and interesting to reformulate our arguments 
in covariant gauges. In addition, the connection of our formalism to the effective action approach 
to small-x and the Regge limit, Refs. [23123 should provide further insight. 
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A Power counting with contracted vertices 

In this appendix we will include the possibility of contracted vertices in the reduced diagram in 
Fig. These are associated with internal lines (collapsed to a point) which are off-shell by ^/s. 
Our analysis closely follows 27j and 

If we go back to the argument that led us to Eq. (fTH|) for the superficial degree of IR divergence 
for the soft part, we see that the same reasoning as in the case of elementary vertices applies to 
the case of contracted vertices since the result ()15|) has been obtained by means of dimensional 
counting. 

The analysis of contracted vertices connecting jet lines only is, however, more subtle. We 
have to demonstrate that the suppression factors corresponding to the contracted vertices are at 
least as great as the ones for the elementary vertices. The expression (j^^ tells us that we can 
restrict ourselves to the two and three point vertices. For these cases, we analyze the full two and 
three-point subdiagrams, by studying the tensor structures that are found after integration over 
their internal loop momenta. 

Before we discuss all the possible structures, we state some results which will be essential for 
the upcoming analysis. The first one is the simple Dirac matrix identity 

^ ^ 2(a-6) ^. (106) 

The other two follow from Eqs. ((Zj) and (fTT|l for the gluon propagator in Coulomb gauge, and 
hold for any jet momenta scaling as ~ ~ -^"j A^/^) collinear to the momentum pA 
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defined in Eq. (0) 



l'^N^p{lA,v) = 0{\^/^^s), (107) 

for all components of (3. We now proceed to discuss the particular cases. 

Ghost self-energy: The most general covariant structure is, using p ■ p = p"^, ^ 

u{p,p) =p■p/(///i^ pVfi^ (108) 

where is a scale introduced by a UV/IR regularization of Feynman diagrams and p is the 
momentum of an internal jet line. Strictly speaking, the covariants should be formed from the 
vectors p and rj, but since p has nonzero light-cone components, we can use Eq. (jHJ, to express rj in 
terms of p. The maximum degree of divergence for the ghost self-energy occurs when the internal 
lines become either parallel to the external momentum p or soft. The most general pinch singular 
surface consists of a subdiagram of collinear lines moving in a direction of the external ghost. 
This subdiagram can interact with itself through the exchange of soft quanta. Power counting 
arguments similar to the ones given in Sec. 13.21 show, however, that there is no IR divergence 
for these pinch singular points. This shows that the dimensionless function / in Eq. p08j) is IR 
finite. Hence the combination [tree level ghost propagator] - [ghost self-energy] - [tree level ghost 
propagator], [l/{p ■ p)]Il{p,p) [l/{p-p)], is suppressed at least as much as a single tree level ghost 
propagator, l/{p ■ p). Therefore the contracted two point ghost vertex within a jet subdiagram 
contributes at least the same suppression as a single tree level ghost propagator. 

Gluon self-energy: With external momentum p, its most general tensor decomposition has the 
form 

^^lv (P, P) = g^Jiu fi + P^,Pu /2 + P^,Pu /s + (Pf^Pu + Pf,Pu) fi ■ (109) 

As verified by explicit one-loop calculations in Refs. P?| and the gluon self-energy in Coulomb 
gauge is not transverse. In Eq. ()109p . the fi = fi (p^//x^, p^//^^, «s(/^)) are dimensionless functions. 
Contracting U^^ with tree level gluon propagators, and using Eq. pO|) . the last two terms in Eq. 
fll09j) drop out and the first and the second terms give at least one factor of in the numerator, 
which cancels one of the (1/p^) denominator factors. Since the maximum degree of IR divergence 
for the gluon self-energy occurs when all the internal lines become either collinear to the external 
momentum p or soft, we can use the results of the power counting of Sec. l3.2l to demonstrate that 
the dimensionless functions fi are at worst logarithmically divergent. Therefore the combination: 
gluon jet line - 2 point gluon contracted vertex - gluon jet line, behaves the same way as a gluon 
jet line for the purpose of the jet power counting. 

Fermion self-energy: In the massless fermion limit, the most general matrix structure of the 
fermion self-energy is 

^iP,p) =f^9i+ P92, (110) 

with dimensionless functions gi = 5'i(p^//U^, as(/i)), i = 1,2. When we sandwich the 

fermion self-energy between the tree level fermion denominators, the first term in Eq. piOp 

^In the rest of this subsection we are concerned the momentum factors only, and we omit dependence on the 
color structure. 
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behaves the same way as the tree level fermion propagator, modulo logarithmic enhancements due 
to the function gi. The second term, however, is absent from the fermion self-energy as was shown 
in Ref. [21] using the method of induction and Ward identities. The idea was to study a variation 
of the fermion self-energy by making an infinitesimal Lorentz boost on the external momentum. 
This implies a relationship between the (r + 1) and the r-loop self energy. Assuming that the term 
proportional to p is absent from the r-loop expansion Sen shows that it is also absent from the 
(r + l)-loop expansion. So the first term in Eq. (jll(J|) is the only possible structure of the fermion 
self-energy when its external momentum is jet like and approaches mass shell. 

Now let us investigate the 3 point functions. 
Fermion- gluon- fermion vertex function: T^, can depend on vectors that scale as Ia, I'a iii Eq. 
(|1(J7|) . provided all momenta external to the contracted vertex are collinear to momentum 
given in Eq. It has one Lorentz index, /i, and neglecting the fermion masses, it contains an 
odd number of gamma matrices. This implies that the most general tensor and gamma matrix 
expansion of involves 

2- 7^t ^ / {Ia ■ Ia) and all permutations of 7^, ^, ^, 

3. hl^A/lA. hll/{lA-lA), VAl^A/ilA-lA), Mil l\- 

The differences between the listed set of structures and other possible combinations are C(A^/^ \fs)i 
as can be shown using Eqs. p()(jj) - p()7j) . The listed gamma matrix structures are multiplied by 
dimensionless functions, which can depend on the combinations l\, f^, l'^, l'^, besides the 
renormalization scale and the running coupling. Using the arguments similar to the ones leading 
to Eq. (j23|) . we easily verify that the above mentioned dimensionless functions are at most 
logarithmically divergent. Next we analyze the possible Dirac structures. 

1. The first case has the same structure as the elementary vertex, and therefore causes the 
same suppression as the elementary vertex. 

2. The fermion-gluon-fermion composite 3-point vertex is sandwiched between the factors 
and I/ai originating from the numerators of the fermion propagators external to the composite 
vertex. Therefore the terms from case 2 where ^ is on the first or third position in the string 
of the gamma matrices provide a suppression On the other hand in the case, when 
^A is in the middle of this string of three gamma matrices, we encounter the combination 
VAlfj. Va after taking into account the numerators of the external fermions. Using Eq. (jlU6p . 
we can immediately recognize that this combination provides a suppression A^^^. 

3. Based on the preceding arguments it is obvious that also the structures included in item 3 
supply at least the same suppression factor as the elementary vertex. 

Therefore, we conclude that the composite 3-point fermion-gluon-fermion vertex behaves as the 
elementary vertex for the purposes of the jet power counting. 

Three gluon vertex: V^,yp, with external momenta collinear to momentum pa- This vertex 
can depend on momenta Ia-, I a defined above and the metric tensor (^q,^. Taking into account the 
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dimension of the 3 gluon Green function, its only possible tensor structure involves combinations of 
the form [fif^^, /^ + perm. + C(A^/^ ^/s)] and [I'X I a ll/^A + ^i^^^'^ v^)]. with all possible replacements 
of Ia^ Ia- These tensor structures are multiplied by dimensionless functions. The former is the 
same as in the case of an elementary vertex and it therefore supplies the same suppression factor 
as the elementary vertex. The latter also provides the same suppression as the elementary vertex, 
since the two momenta, say after being contracted with the propagators of the external 

gluons, give suppression factors, as in Eq. ()l()7j) . which cancel the 1//^ enhancement. The leftover 
momentum provides the same suppression factor as the elementary vertex. Using the collinear 
power counting of Sec. 13.21 one can immediately see that the IR divergence of the dimensionless 
functions multiplying these tensor structures is not worse than logarithmic. Hence, there is a 
suppression factor A^/^ associated with every contracted 3 gluon vertex. 

Ghost - gluon - antighost three point vertex: When all lines external to the contracted vertex 
are of the order I a, the most general tensor structure for this contracted vertex is 

/^/ii + p;/i2 + o(A'/'v^), (111) 

with dimensionless functions hi = hi (/^//i^, T^Z/i^, as{^)), i = 1,2, which are at most logarithmi- 
cally IR divergent. Using Eq. (|1U7|) . we see that when the momenta in Eq. are contracted 
with the tree level gluon propagator, we get a suppression of the order of the transverse jet mo- 
mentum, and that this contracted vertex gives the same suppression as the elementary three point 
vertex, at least. 



B Varying the Gauge-Fixing Vector 

In this appendix we study the effect of an infinitesimal boost, performed on the gauge fixing vector 
rj, on an expectation of a time ordered product of fields, denoted by O, taken between physical 
states. The gauge-fixing and the ghost terms in the QCD lagrangian are 

^g.iXx) = -^gl{x), 

jCghost{x) = -ba{x)6BRSga{x) /6A, (112) 

respectively. In Eq. flll2j) . 6 A is a Grassmann parameter defining the BRS transformation, ba{x) 
is an antighost field and 

gaix) = -d ■ Aaix) = -id-ir]-d)ri)- A„(x). (113) 

Let us consider an infinitesimal boost with velocity 6(3 on a gauge fixing vector rj performed in 
the plus-minus plane 

T] ^ T]' = r] + fi6p, (114) 

where the vectors rj and rj are defined in Eqs. Q and (jSHI), respectively. Since only the gauge 
fixing and the ghost terms in the QCD lagrangian depend on rj, we can write to accuracy (9(5/5^) 

6 <0> = < 0{r]') > - < 0{r]) > = < r^^^^ 6(3 > 
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d^x < 0{r]) ga{x) 6ga{x) > -i j d^x < 0{r]) ba{x) 6 {6brs 9a{x)/6A) > . 

(115) 



Using the BRS invariance of the QCD lagrangian and the BRS transformation rule for an antighost 
field 

5BKsba{x)/5A = ^ga{x), (116) 



we arrive at 



5<0>=-i j d^xK (5brsO/5A) ha{x) 5ga{x) > . (117) 
Taking a variation of ga{x) in Eq. ()117|1 . we obtain 

5<0>=-t [d^x< (5brs 0/6A) baix) {{fi ■d)7] + {r]-d)fj)- Aa{x) > . (118) 



Substituting for O a product of n gluon fields, we can use Eq. ()118|) . together with the rule for 
the BRS transformation of a gluon field 

Sbrs Alix)/5A = d.c'^ix) + g^r^'^Alixyix), (119) 

with c°'{x) representing the ghost field, to derive the gauge variation for a connected Green func- 
tion. However, our jet functions are one-particle irreducible in external soft lines and we therefore 
cannot apply Eq. (jllSp directly, and must find an analog for this subset of diagrams. The modi- 
fication of Eq. ()118|) due to the restriction to IPI diagrams is, however, not difficult to identify. 

Let us consider the graphical analog of the derivation of Eq. ()118j) just outlined. The variation 
in rj may be implemented as a change in the gluon propagator and, in Coulomb gauge, the ghost- 
gluon interaction, which is also ?7-dependent. This is the viewpoint that was taken in axial gauge 
in Ref. P2j- At lowest order in the variation, the modified gluon propagator produces scalar- 
polarized gluon lines, which decouple through repeated applications of tree-level Ward identities 
to the sum over all diagrams. The relevant tree-level identities are given in p3|. We need not 
describe these identities in detail here. We need only note that they are to be applied to any 
diagram in which a scalar polarized gluon appears at an internal vertex. Every such application 
produces a sum of diagrams, each of which fall into one of two sets: 1) diagrams in which an 
internal gluon line is transformed to a yet another ghost line ending in a scalar polarization, and 
2) diagrams in which one gluon line is contracted to a point. The new vertex formed in the former 
case is the ghost term, and in the latter case it is the ghost-gluon vertex of the BRS variation 
()119j) . Eq. ()118|1 must result from the cancellation of all diagrams, set 2), in which an internal 
gluon line is contracted. Contracted external lines provide the ghost-gluon terms, and the ghost 
lines of set 1) eventually provide the ghost terms of the BRS variations ()119|) of external fields in 
Eq. (fTT8|) . 

The simplicity of the tree level Ward identities puts strong limitations on the sets of diagrams 
that can combine to form different diagrammatic contributions to Eq. pi8j) . For diagrams of 
set 1), the topology of the original diagram is unchanged, and a IPI diagram remains IPI. For 
diagrams of set 2), generally IPI diagrams remain IPI, except in the special case of a diagram that 
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a) b) c) 

Figure 12: Two-loop diagram illustrating the idea of the tulip-garden formalism. Ti,T2,T3 are 
the possible tulips. 

is two-particle reducible, with these two lines separated by a single propagator. In this case, the 
contraction of the internal line that separates the other two will bring those two lines together at a 
single vertex, producing a diagram precisely of the topology shown in Fig. HJ On the one hand, by 
Eq. ()118|) all such diagrams must cancel in the full perturbative sum. On the other hand, the same 
topology results from a diagram that is one-particle reducible with respect to a single line, which 
is then contracted as a result of the tree-level Ward identity. The latter diagram, however, is not 
included in the set of IPI diagrams with which we work. The application of the Ward identity 
to IPI diagrams only, therefore, results in terms that would cancel this special set of one-particle 
reducible diagrams, in which the only line that spoils irreducibility is contracted to a point. These 
are the diagrams shown in Fig. in which the ghost-gluon vertex of Eq. ()119p is inserted between 
one-particle irreducible subdiagrams in all possible ways. The ghost line ending at this composite 
vertex is continuously connected to the variation of a gluon propagator, according to Eq. ()118|) . 
The full composite vertex of the Ward identity in Eq. (jll8|) appears only at true external lines of 
the IPI jet. This vertex is given by the momentum factor in Eq. (|37p and is represented by the 
double line crossing a gluon line in Fig. IT^ below. Diagrams that are reducible in one or more 
internal lines can be treated in a similar manner. The "left-over" terms in the Ward identities 
for each set of diagrams of definite reducibility properties (IPI, 2PI, etc.), must cancel in the full 
sum, reproducing the identity for Green functions, Eq. (jll8|) . 

C Tulip- Garden Formalism 

In this appendix we illustrate how a given Feynman diagram contributing to the process in the 
leading power can be systematically written in the form ()27j) . For concreteness let us consider a 
two loop diagram where the quarks interact via the exchange of a one rung gluon ladder as in Fig. 
IT^ The important contributions of this diagram come from the regions when all of the exchanged 
gluons are soft. Fig. IT^ or when the gluons attached to the A quark line are soft, while the rest 
of the gluons carries momenta parallel to the — direction (they belong to jet B), Fig. IT^ . or 
when the two gluon lines attached to the B quark line are soft and the other gluons are collinear 
to the -|- direction (they belong to jet A), Fig. IT^ . The possible central soft exchange parts are 
called tulips. In our case the possible tulips are denoted as Ti,T2,Ts in Fig. [T21 The garden is 
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defined as a nested set of tulips {Ti, . . . , T„} sucli tliat Tj C Tj+i for i = 1, . . . , n — 1. In Fig. [T21 
{Ti}, {T2}, {T3}, {Ti.Tg}, {Ts.Tg} are the possible gardens. 

For a given tulip we make the soft approximation, consisting of attaching a soft gluon to jet 
A via the — component of its polarization only and to jet B via the the + component of its 
polarization. The result of this soft approximation for a given Feynman diagram F corresponding 
to a tulip T is denoted S{T)F. It has obviously the form of Eq. (j27|) . Following the prescription 
given in Refs. [22] and |3| we write the contribution to a given diagram F in the form 

F = ^(-l)"+i5(ri) . . . 5(r„) F + Fn, (120) 

G 

where the sum over inequivalent gardens, as defined bellow, G in Eq. (|12Up is understood. The 
meaning of this expression is the following. For a given garden consisting of a set of tulips 
{Ti, . . . ,Tn}, we start with the largest tulip T„ and make the soft approximation for the gluon 
lines coming out of it. Then for T„_i we proceed the same way as for T„. If some of the lines 
coming out of T„_i are identical to the ones coming out of T„ we leave them untouched. For 
instance, if we consider a garden {T2,T^} from Fig. ^[ we first perform the soft approximation 
on tulip T3 and then proceed to tulip T2. However the lines coming out of T2 and T3 which attach 
to the B quark line are identical so when performing S{T2)S{Ts)F we leave these gluon lines out 
of the game and make soft approximations only on the gluon lines attaching to the ladder's rung. 
Two gardens are equivalent if the soft approximation is identical for both of them. is defined 
by Eq. (|12m) . The contribution to F^ comes from the integration region where \k\ > ^/s for 
all gluons coming out of the central soft part. As a result, the contribution to is suppressed 
by positive powers of ^/—i/^/s. Therefore we can ignore the contribution from Fji within the 
accuracy at which we are working. 
We can now rewrite Eq. ()120|) . as 

T Vg,t„=t / 

This expression is indeed in the form of Eq. ()27j] since the term S{T)F is of that form and the 
subtractions ^(— l)""'"^5'(Ti) . . . S(Tn-i) modify only the soft function S in Eq. (j2Z|), but do 
not alter the form of the equation. We can therefore conclude that the contribution to a given 
Feynman diagram in leading power can be expressed in the first factorized form given by Eq. ()27|) . 

D Feynman Rules 

In Fig. ^1 we list the Feynman rules for the lines and the vertices encountered in the text. The 
double lines are eikonal lines, while the dashed lines represent ghosts. The four vectors 1], fj are 
defined in Eqs. (jHI) and respectively. The conventions for the gluon-ghost and gluon-eikonal 
vertices (third and second from the bottom of Fig. IT!?|) are the following. We start with a color 
index of a gluon external to the diagram defining the evolution kernel, see for instance Fig. |H^, 
then proceed to the gluon internal to the diagram and finally to the ghost /eikonal line in order to 
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Figure 13: Feynman rules for the eikonal lines, ghost lines and special vertices. 

assign the color indices of fate- For the three point antighost - gluon - ghost vertex at the bottom 
of Fig. ^[ we start with an antighost (arrow flowing out of the vertex) then proceed to the ghost 
and finally we reach the gluon line. 

E Origin of Glauber Region 

In this appendix we exhibit the origin of the Glauber (Coulomb) region using the two-loop diagram 
shown in Fig. El Consider a situation when the upper gluon loop is a part of J a- Momentum k 
of the exchanged gluon flows through jet lines with momenta I2 = I — k and = pa — I — q + k. 
The components of k can be pinched by double poles coming from the denominators of the gluon 
propagators k"^ + ie and (g — k)"^ + ie. In addition to these pinches, the component k~ can be 
pinched by the singularities of the jet lines I2 and I3, at values 

k- = r-^^, 

k- = + (122) 
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Figure 14: Two loop diagram demonstrating the origin of the Glauber (Coulomb) region. 

The two poles given by Eq. ()122|) are located in opposite half planes since in the region considered 
I2, It > 0. This indicates that we must consider the possibihty that the different components 
of the soft momentum k can scale differently. For instance, we can have fc"*" ~ fc^ ~ o^^/s and 
k~ ~ Xy/s where A ^ a ^ 1. Indeed, the power counting performed in Sec. 13.21 shows that 
the singularities originating from these regions can produce a logarithmic enhancement. We also 
note that it is only minus components that are pinched in this way by the lines in J a, and plus 
components by the lines in Jb- 
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